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We investigate nonequilibrium properties of the single impurity Anderson model by means of the functional 
renormalization group (fRG) within Keldysh formalism. We present how the level broadening T/2 can be used 
as flow parameter for the fRG. This choice preserves important aspects of the Fermi liquid behaviour that the 
model exhibits in case of particle-hole symmetry. An approximation scheme for the Keldysh fRG is developed 
which accounts for the frequency dependence of the two-particle vertex in a way similar but not equivalent to a 
recently published approximation to the equilibrium Matsubara fRG. Our method turns out to be a flexible tool 
for the study of weak to intermediate on-site interactions U < 3r. In equilibrium we find excellent agreement 
with NRG results for the linear conductance at finite gate voltage, magnetic field, and temperature. In nonequi- 
librium, our results for the current agree well with TD-DMRG. For the nonlinear conductance as function of 
the bias voltage, we propose reliable results at finite magnetic field and finite temperature. Furthermore, we 
demonstrate the exponentially small scale of the Kondo temperature to appear in the second order derivative of 
the self-energy. We show that the approximation is, however, not able to reproduce the scaling of the effective 
mass at large interactions. 
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I. INTRODUCTION 

Due to the enormous experimental progress in the inves- 
tigation of nanoelectronic and molecular systems, one of 
the central issues of theoretical condensed matter physics 
is the understanding of nonequilibrium phenomena in small 
strongly correlated quantum systems coupled to an environ- 
ment. Traditionally such systems were investigated in the 
context of dissipative quantum mechanics, where energy ex- 
change of a few-level system with a single bath of phonons 
was considered. 1 With the realization of various devices in- 
volving metallic islands 2 , quantum dots 3 , single molecules 4 , 
or quantum wires*^, a tunneling coupling to leads was real- 
ized, which enabled the controlled study of local quantum sys- 
tems coupled via particle exchange to the environment. Fur- 
thermore, in the presence of several leads, it became possible 
to study such systems in the presence of a finite bias voltage or 
temperature gradients (inhomogeneous boundary conditions), 
i.e. a nonequilibrium and current-carrying state can be real- 
ized in the stationary situation. 

Whereas in equilibrium there are powerful numerical and 
analytical techniques to describe linear transport or spectral 
properties of mesoscopic systems, the available methods in 
nonequilibrium are quite restrictive and still under develop- 
ment in the nonperturbative regime. Exact solutions exist 
for some special cases 7 " 9 . Scattering Bethe-Ansatz solutions 
have been applied to the interacting resonant level model 10 
and the Anderson impurity modelii, but the full understand- 
ing of the stationary state in all regimes has not yet been 
reached. To describe the steady state in the presence of a finite 
bias voltage, numerical renormalization group (NRG) meth- 
ods using scattering wavesi^, time-dependent density ma- 
trix renormalization group (TD-DMRG) 13 , quantum Monte 
Carlo (QMC) with complex chemical potentials 14 , and an ap- 
proach based on iterative summation of path-integrals (ISPI) 



(Ref. [I5I) have been developed, but the efficiency of these 
methods is still not satisfactory in the regimes of strong 
Coulomb interaction and/or large bias voltage. To calculate 
the time evolution into the stationary state, numerical tech- 
niques such as TD-NRGr&i^ TD-DMRG^i^ ISPI^ and 
QMC in nonequilibrium 19,20 have been used, as well as a Fock 
space formulation of the multilayer multiconfiguration time- 
dependent Hartree theory 2 ^. However, the description of finite 
bias or the long-time limit still remains difficult. 

Perturbation expansions in the coupling between the local 
system and the reservoirs or in the Coulomb interaction on 
the local system can not cover the most interesting regime of 
quantum fluctuations and strong correlations at very low tem- 
peratures. To improve them significantly, the most promising 
analytic approaches are perturbative renormalization group 
(RG) methods for nonequilibrium systems. Several meth- 
ods have been proposed, which differ significantly concern- 
ing the extent up to which nonequilibrium aspects are taken 
into account and how the perturbative expansion in the renor- 
malized vertices is set up. One of the first formally exact 
nonequilibrium RG methods for zero-dimensional quantum 
systems (quantum dots) was developed in Refs. |22|424| and is 
called the real-time renormalization group method (RTRG). 
This technique is based on an expansion in the renormalized 
coupling between quantum system and reservoir, whereas the 
correlations on the local system are taken exactly into ac- 
count. Later on, this method was formulated in pure fre- 
quency space and combined with a cutoff procedure on the 
imaginary frequency axis (introduced in Ref. l25h . the so- 
called RTRG-FS method 26 . Within this technique it was pos- 
sible to solve analytically the generic problem how RG flows 
are cut off by the physics of transport rates. Furthermore, 
it was demonstrated how all static and dynamic properties 
can be calculated analytically for the nonequilibrium Kondo 
model in weak coupling 27 " 29 and for the interacting reso- 
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nant level model (IRLM) in the scaling limit , which are 
two fundamental models for the physics of spin and charge 
fluctuations, respectively. Another perturbative RG method 
in nonequilibrium (called PRG-NE), which also expands in 
the reservoir-system coupling, was developed in Refs. l3lll32l 
where the slave particle approach was used in connection with 
Keldysh formalism and quantum Boltzmann equations. In 
these works it was investigated for the first time how the volt- 
age and the magnetic field cuts off the RG flow for the Kondo 
model and how the frequency dependence of the vertices influ- 
ences various logarithmic contributions for the susceptibility 
and the nonlinear conductance. A real-frequency cutoff was 
used and the RG was formulated purely on one part of the 
Keldysh-contour disregarding diagrams connecting the upper 
with the lower branch. This procedure turns out to be suffi- 
cient for the Kondo model to calculate logarithmic terms in 
leading order but the cutoff by relaxation and decoherence 
rates was included intuitively (recently an improved version 
of this method included parts of the spin relaxation and deco- 
herence rates via self-energy insertions 33 ). An alternative mi- 
croscopic approach to RTRG-FS for combining relaxation and 
decoherence rates within a nonequilibrium RG method for the 
Kondo model was proposed in Ref. 34, where flow-equation 
methods 35 were generalized to the nonequilibrium situation. 
Within this method, it was shown for the Kondo model that 
the cutoff of the RG flow by spin relaxation/decoherence rates 
occurs due to a competition of 1-loop and 2-loop terms on the 
r.h.s. of the RG equation for the vertex. 

Whereas all of these perturbative RG methods are very suc- 
cessful, they have two important drawbacks. First, they ex- 
pand all in the reservoir-system coupling and thus can only be 
used reliably in the regime of weak spin or orbital fluctuations 
(strong charge fluctuations seem to be covered by RTRG-FS 
as was demonstrated recently for the IRLM in Ref.[30h. Sec- 
ondly, the RTRG-FS and the PRG-NE scheme work in a basis 
of the many-particle eigenfunctions of the isolated local quan- 
tum system. As a consequence larger systems such as multi- 
level quantum dots, large molecules, or the crossover to the 
quantum wire case cannot be addressed due to an exponen- 
tially large number of relevant many-particle states. There- 
fore, another class of perturbative RG methods in nonequi- 
librium have been developed which are based on an expan- 
sion in the renormalized Coulomb interaction on the local 
quantum system by combining the Keldysh formalism with a 
quantum-field theoretical formulation of functional RG within 
the 1-particle irreducible Wetterich scheme 36 . In the follow- 
ing we call this approach the fRG-NE method. It has been 
proposed and applied to transport through quantum dots and 
quantum wires in Refs. Similar implementa- 

tions of the method have been developed for the application 
to bulk systems, investigating for instance quantum critical- 
ity in nonequilibrium 40 or the long-time evolution of a Bose 
gas starting from a nonequilibrium state 41 . Since the method 
is non-perturbative in the reservoir-system coupling and is 
parametrized in terms of single-particle states, it has the ad- 
vantage that the case of low Voltage and Temperature can be 
described and it can be applied to multi-level quantum dots 
and quantum wires, at least if Coulomb interactions are mod- 



erate. 

So far, the fRG-NE scheme was used in approximation 
schemes with frequency independent self-energy. This ap- 
proach describes reliably the leading effects of weak interac- 
tions in an effective single-particle picture. At higher interac- 
tions the decay rates generated by inelastic effects of the inter- 
action are no longer negligible and may drastically change the 
behaviour of the system. These effects can only be described 
on the level of a frequency dependent self-energy, which is 
the subject of the present paper. We will address the problem 
of enhancing the level of the truncation scheme used in fRG- 
NE, in order to enable it to describe interaction induced decay 
rates. The flow equations for the vertex functions are of such a 
structure that a frequency dependent self-energy can only re- 
sult from a frequency dependent two-particle vertex. This is in 
turn a very complicated object, depending on four frequency, 
state and Keldysh indices. Conservation laws for spin, mo- 
mentum, and frequency, and the symmetry properties of the 
vertex function discussed in Ref.|42jcan be used to reduce the 
complexity of the problem. Nevertheless, the two-particle ver- 
tex function remains so complicated that it seems preferable 
to develop the methodology first for a physical system with as 
few interacting degrees of freedom as possible. This restricts 
at least the dimensionality of the state dependence of the ver- 
tex function. Insight gained at this level can later be used to 
tackle larger systems. 

We choose to investigate the single impurity Anderson 
model 43 (SIAM) which is considered to be a generic model 
for strong local correlations 44 and is a minimal model to study 
the interplay of charge and spin fluctuations. It consists of 
a single-level quantum dot with two spin states and on-site 
Coulomb repulsion U, which is coupled to leads. Experimen- 
tally, it can be realized for small quantum dots or molecules. 
In the Coulomb-blockade regime, this model is equivalent to 
the Kondo model 45 and the unitarity limit of the Kondo effect 
has been measured 46 . The perturbation theory of this model 
in U is completely regular 47 " 50 . However, the low energy be- 
haviour at large interaction U is governed by the so called 
Kondo scale 44 which is exponentially small in U and hence 
cannot be described in any finite order perturbation theory. 
This makes the model an interesting object for RG studies. 
Since the model can be described in equilibrium very accu- 
rately by the numerical renormalization group (NRG) 5 -^ and 
even exact results for thermodynamic properties are available 
from Bethe-Ansatz 52 there is the possibility to benchmark our 
results at least in equilibrium. 

In case of thermal equilibrium the Matsubara fRG has so far 
been applied to the SIAM in two distinct formulations. One of 
them is a recent study based on Hubbard-Stratonovich fields 
representing spin fluctuations 53 . The fRG-NE implementation 
presented in this paper does not pursue this idea but is based 
on a weak coupling expansion analogously to the second class 
of equilibrium fRG studies of the SIAM 54 " 57 . These studies in 
turn have been done in two different truncation schemes. One 
of those schemes reduces the flow of the two-particle vertex 
to a renormalization of the static interaction strength and pro- 
duces a frequency independent self-energy. This approxima- 
tion is able to reproduce a Kondo scale exponentially small in 
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U which appears e.g. in the pinning of the level to the chem- 
ical potentia l 54 ' 55 . It turned out to be a quite powerful tool 
for the description of static properties of diverse quantum dot 
geometries even at fairly large interactions 54 . The approxima- 
tion is however unable to describe finite frequency properties. 
Its restriction becomes manifest very clearly in the shape of 
the spectral function: the approximation of a static self-energy 
leads to a Lorentzian resonance peak. The sharp Kondo reso- 
nance, the side bands, the suppression of the Kondo resonance 
with temperature are features which cannot be described in 
principle by the static approximation. 

The second more elaborate approximation scheme which 
has been used within the weak-coupling equilibrium fRG 
keeps the frequency dependence of the two-particle vertex 
function and neglects the three-particle verte x 56 ' 57 . It yields 
quantitatively good results for small to intermediate interac- 
tion strengths. While finite frequencies properties are acces- 
sible to this frequency dependent fRG approximation for not 
too large interactions, it has been observed that the large inter- 
action asymptotics of this refined version are worse than the 
ones found in the static fRG; no Kondo scale exponentially 
small in U has been found 57 . Also a technical disadvantage 
related to the use of the Matsubara formalism becomes ap- 
parent: while the self-energy data are quite accurate on the 
imaginary frequency axis for small and intermediate U, the 
analytic continuation to the real axis was instable for nonzero 
temperatures, working only for certain parameter sets. There- 
fore only observables which are accessible from the imaginary 
frequency axis (without resorting to analytic continuation) can 
be systematically studied at finite temperature 57 . This find- 
ing is an additional motivation to study a frequency dependent 
fRG within the Keldysh formalism, being formulated on the 
real frequency axis from the outset. 

Apart from the possibility to access real frequency proper- 
ties at finite temperature we envisage the opportunity to an- 
alyze the non-equilibrium behaviour of the SIAM. This has 
been recently in the focus of diverse publications. Among 
them is also a non-equilibrium fRG study that is based on a 
real-frequency cut-off and a static approximation scheme 38 ' 39 ; 
results known from the Matsubara fRG could be partially re- 
produced by this method. It suffers however from the viola- 
tion of causality related to the choice of the real-frequency cut- 
off. References 581:59 present perturbative studies that apply 
to moderate interactions and the special situation of particle- 
hole symmetry and vanishing magnetic field. The RTRG 
method was used in Ref.|22jto describe the mixed valence and 
empty-orbital regime of the SIAM at finite bias. The results 
seem to be reliable but it was not possible to study the Kondo 
regime since essential processes describing spin fluctuations 
were neglected. In Ref. [l4j a time independent description 
of non-equilibrium based on a density operator for the steady 
state is used to make the problem accessible to QMC. The 
main challenge of this method is the numerical analytic con- 
tinuation of two imaginary quantities to the real axis. A recent 
improvement of nonequilibrium QMC applied to the SIAM at 
zero magnetic field and zero gate voltage 20 has given reliable 
results for U < 3 — 5F = 6 — 10A (in our notation A = T/2 
denotes the level broadening, i.e. T is the full width at half 



maximum of the noninteracting spectral density) for the dot 
occupation and the nonlinear current. The ISPI approach 15 
provides an access to the current through the system for mod- 
erate U and not too low temperatures. A very promising tool 
is the recently introduced scattering states NRG 12 . Unfortu- 
nately the results for the conductance obtained by this method 
still bear a considerable numerical uncertainty. Spataru et al. 60 
generalized the GW approximation to non-equilibrium and 
concentrate on the Coulomb blockade regime of the model. 
The TD-DMRG technique^ tries to access the steady state 
transport features of the model from the transient regime. Data 
of the ISPI and the TD-DMRG methods for the current at 
moderate interactions have been found to agree very well with 
our fRG results.^ 

In contrast to many other methods, we propose in this pa- 
per a very flexible approach which provides reliable results 
for moderate interactions U <3T ~ 6A. In equilibrium, we 
will show that the linear conductance agrees very well with 
NRG data at finite gate voltage, magnetic field, and tempera- 
ture. In nonequilibrium at finite bias V, the nonlinear current 
I(V) is found to agree very well with TD-DMRG data 13 . This 
provides the main evidence that our results can be trusted also 
for the nonlinear conductance G(V) at finite magnetic field B 
and finite temperature T. In particular, at T = B = 0, we find 
that our results for G(V) do not show anomalous peaks as ob- 
tained within fourth order perturbation theory 59 . Furthermore, 
we show that even the exponentially small scale of the Kondo 
temperature can be identified in the second order derivative of 
the self-energy and the various Fermi-liquid relations are re- 
produced. However, the effective mass still does not contain 
an exponentially small scale leading to a too large broadening 
of the spectral density for interaction strengths U > 3T, in line 
with the equilibrium Matsubara fRG 5 - 7 -. 

The paper is organized as follows. SectionHJintroduces the 
model and its treatment within Keldysh formalism. In sec- 
tion [Hi] we recapitulate the core elements of the fRG for irre- 
ducible vertex functions. A motivation and discussion of the 
choice of hybridization as flow parameter follow in sectionlTVl 
In section [V] we describe the more basic frequency indepen- 
dent approximation of the flow equations and show that in 
the case of equilibrium and zero temperature we reproduce 
exactly the flow equations known from the Matsubara fRG 
Section|VT]is then devoted to the fRG in frequency dependent 
approximation. There are three channels contributing to the 
flow of the two-particle vertex, and all three need to be taken 
into account. We describe an approximation scheme that sim- 
plifies the functional form in which the two-particle vertex 
depends on frequencies. In section [VTTl we demonstrate that 
the fRG approach in the chosen truncation and approximation 
scheme preserves the Fermi-liquid relations for the (imagi- 
nary part of the) self-energy; furthermore we determine the 
fRG estimate for the Fermi-liquid coefficients. Section IVIIII 
finally presents numerical results obtained from the fRG and 
compares them to other methods. A conclusion is given in sec- 
tion [DC] The Appendices present some mainly technical con- 
siderations. In Appendix |A] we identify a set of independent 
components which completely determine the two-particle ver- 
tex function. The precise form of the flow equations in terms 
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of these components is given in Appendices iBl and ICl In Ap- 
pendix[D]we recover within the fRG approach the nonequilib- 
rium Fermi-liquid relation known from Ref. |62j. 



H. SIAM AND KELDYSH FORMALISM 

The single impurity Anderson model 43 under consideration 
consists of a single electronic level with on-site repulsion, 
which is coupled to two noninteracting reservoirs addressed 
as left (L) and right (R). We denote the single particle states 
of the impurity by a =^,|= according to the state 

of the electron spin. The matrix element of the on-site two- 
particle interaction is 



{o[02\v\aio 2 } 



U, if o[ = <7i 
0, else, 



o 2 , 



(1) 



with U > 0, where we used the notation a = —a. In standard 
notation of second quantization the Hamiltonian is given by 



H — //dot + //res + //( 



coup? 

u 



(2a) 



//dot = £ [eV g - oB - — ) d\d c + Udld^dU^, (2b) 



=L.R 



(2c) 



//coup = £#ioup - L L / d ^ (Vfe4c*.a + h-c), (2d) 



where the annihilators t/ a ,Q r o- an d creators 4,ct „ obey the 
usual fermionic anti-commutation rules. The single-particle 
energies of the dot 



£a = eV g -oB-U/2 



(3) 



depend on the gate voltage V s and the magnetic field B and are 
shifted by (-U/2) such that particle-hole symmetry is given 
when eVg equals the chemical potential. 

We are interested in the stationary state which emerges a 
long time after the system has been prepared in a product den- 
sity matrix 



P (t ) = Po = PL ® Pdot ® PR, fo 



(4) 



Here, Pl and Pr describe each an individual grand-canonical 
equilibrium characterized by the Fermi function 



fr((0) 



1 



L,R. 



(5) 



(We use units with h = 1 and &b = 1 throughout this paper.) 
The temperature T is assumed to be equal in both reservoirs, 
whereas a possible difference between the two chemical po- 
tentials accounts for a finite bias voltage, 



eV = /i L -iJ, R . 



(6) 



We measure single-particle energies relative to the mean 
chemical potential by setting 



(7) 



We describe the system in the framework of Keldysh 
formalism 63 " 66 , where the single-particle propagator between 
two states q' and q has four components 

G-}-(t\t') = G c qW (t\t') = -iTrp JS('K'( ? ') < 8a > 

(8b) 
(8c) 



G ,|7( f 10 = G «w( { V) = iKpoal(t')a q (t) 



GffiW) = G ^/( f l0 = -ttrpoa q (t)al(t') 



G+ awW) = Gl w W) = -iTrp ^(0flI,(f') (8d) 



called chronologic, lesser, greater, and anti-chronologic, re- 
spectively. Here ST (,^) denotes the time (anti-time) ordering 
operator and a q (t) — d a (t) or Q r(T (f) is in the Heisenberg- 
picture at time t. In the time translational invariant stationary 
state G(t\t') = G(t — f'|0) depends only on the difference of 
the two time arguments and we use the Fourier transform 



G J V(co)= dte^G^m, j,f = T . 



(9) 



Instead of the contour basis with indices j = =p we use the 
Keldysh basis 67 with indices a = 1,2. The transformation to 
this basis is given by 



(10) 



where 



D -n =D M2 = {D -y\- = (D -i )2 | T = J_ (lla) 



D 



+ i 



(D~ ) 



i\i + 



1 

V2- 



(lib) 



The resulting components of the single particle Green func- 
tion 

G^=0, G 1|2 = G Av , G 2|1 =G Ret , G 2|2 = G K (12) 

are called advanced, retarded, and Keldysh, respectively. The 
self-energy is transformed correspondingly, 



(13) 



/ J==F 



leading to 



E 1|1 =E K E l|2 = E Ret ; E 2|1 =E AV ; E 2|2 _ q (14) 

The influence of the reservoirs on the dot Green function 
can be described by reservoir self-energy contributions 68 



l£ a ' (a)) = lAk r V kr gl^((0)VZ 



(15) 



with I = 2, 2 = 1, and g^. being the free propagator in the 
reservoir state k r which is assumed to be spin-independent. 
Explicitly, the individual Keldysh components are 



E&a >) = — I dm' , . 

2% J ca-ca' + iri 

l r -\Av / r \Ret 



Ret 



1 r. , T r (af) 



(16a) 



(16b) 



nZUa) = -i[i -2/ r ((»)]r r ((») ) (i6c) 



5 



T,Hl 




eV 



FIG. 1: Sketch of the model with single particle energy in vertical 
direction. 



where we made use of the hybridization function 

T r {co) = 2n Jdk r \V kr \ 2 d(co-e kr ). (17) 

As we are neither interested in the influence of the band struc- 
ture of the reservoirs nor in the momentum dependence of 
the hopping elements, we linearize the reservoir dispersion, 
% = w r k r (with w r = de kr /dk r ) being the inverse density of 
states at the Fermi level, k r = kf.) and set the hopping to mo- 
mentum independent constants, V kr=k v = V r . As a result, the 
hybridization functions are frequency independent constants, 



r r («) 



2% 



W r 



(18) 



The results presented later apply to the special case of sym- 
metric coupling, Tl = Tr, which is technically easier for the 
implementation of our approximations. 
Defining the total hybridization T via 



r = r. 



(19) 



we find the reservoir-dressed noninteracting dot propagator 
1 



«r(a>) = - .. , ..., v 

CO — £<j + ll /2 

g^(co) = [1 - 2/ eff (to)] [g^(co) -f£(a>)] , 
with the effective distribution function 



CO 



(20a) 

(20b) 
(20c) 



(21) 



A sketch of the model can be found in Fig.Q] 

The current from the left reservoir through the dot to the 
right reservoir is given by§&£2 

/ = «^/dfl)[/L(fl»)-^l(fl))]i;pa(fl)) (22) 



where the spectral density can be obtained from the interacting 
single-particle dot Green function as 



Po(CO) 



1 



-ImG 



Ret, 



CO 



(23) 



In the case Tl = Tr, which we focus on later, the current is an 
odd function of the bias voltage: according to equations (01- 
0, inverting the sign of V means interchanging /l and /r. 
As a consequence, the effective distribution function ( f2TT >. the 
reservoir dressed propagator ( |20l i, the interacting Green func- 
tion, and the spectral function d2"3l remain unchanged, while 
the sign of the current (1221 changes. 

We note that the Hamiltonian is of such a type that the dot 
Green function in the spin basis exhibits the special behav- 
ior under time reversal described in section 5.5 of Ref. l42l 
Therefore the Kubo-Martin-Schwinger conditions which con- 
nect the components of the Green and vertex functions in ther- 
mal equilibrium and which we apply later on take the form of 
a generalized fluctuation dissipation theorem as stated in sec- 
tion 5.6 of that reference. 



III. VERTEX FUNCTIONS AND FUNCTIONAL 
RENORMALIZATION GROUP 

A treatment of the problem within the fRG is set up by mak- 
ing the bare propagator g depend on a flow parameter A . Most 
commonly A is chosen to suppress low energy degrees of free- 
dom. Being functionals of the bare propagator, the interacting 
Green and vertex functions acquire a dependence on A as well, 
which is described by an infinite hierarchy of coupled flow 
equations. 70 The flow parameter is introduced in such a way 
that two values of A are of particular importance: at A = A s tart 
the Green or vertex functions can be determined exactly or 
in reasonable approximation; at A = A s top the free propagator 
takes its original value, g(A s t p) = g, and so do the interacting 
Green or vertex functions. Therefore the interacting Green 
or vertex functions can be found by (approximately) comput- 
ing their flow from A sta rt to A stop . In this paper we focus on the 
flow of the one-particle irreducible vertex functions 36 - 71 which 
has proven to provide a successful approach to the physics of 
diverse low-dimensional correlated electron problem s 72,73 . 

The one-particle irreducible vertex functions can be de- 
rived from generating functionals; details on this approach can 
be found for Matsubara formalism, e.g., in Ref. 74, and for 
Keldysh formalism in Ref. 38. Equivalently, the one-particle 
irreducible n-particle vertex function 



Yin \ W\ i ■ 

'q'y-q'nlqi-qn y l ' 



.,C0'„\C0i,...,C0„) 



(24) 



can be defined diagrammatically as the sum of all one- 
particle irreducible diagrams with n amputated incoming lines 
(having states q\,...q n , Keldysh indices a.\,...,a„ and fre- 
quencies C0i , . . . , C0„) and n amputated outgoing lines (having 



states q\ 



Keldysh indices a[ , . . . , a! n and frequencies 



Co[ , . . . , Co' n ). Due to frequency conservation only (2« — 1 ) of 
the 2n frequencies are independent. By the expression 
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we refer to a function depending only on (2n — 1 ) frequency Here we use shorthand notation like 
arguments, the last one being redundant. 

The vertex functions are anti-symmetric under exchange of 
particles. In fact, the underlying diagrammatics is meant to 
be of the Hugenholtz type 74 , thus based on anti-symmetrized 
interaction vertices 



ii.N^sir^^h^), (29) 



'q l q 2 \q [ q2 




a' 7 is odd, 



(25) 



where q,q2 iq'i^q'j) are incoming (outgoing) single particle 
states and CC\ , OC2 (oc( , OC^) incoming (outgoing) Keldysh in- 
dices. 

In order to evaluate a specific diagram contributing to a ver- 
tex function, one determines the symmetry factor S, the num- 
ber « e q of equivalent lines, the number ni oop of internal loops 
and the permutation P that describes which incoming index is 
connected to which outgoing one. The value of the diagram is 
then given by 



(_l)Kioo P (_i)f nn 

It 



n-l 



n^ln 



(26) 



where one has to sum over all internal state and Keldysh in- 
dices and to integrate over all independent internal frequen- 
cies. Details can be found in Ref . l42l Equation d26l i defines 
the prefactor of the vertex functions in such a way, that they 
can be used themselves as vertices in diagrams with the iden- 
tical prefactor rules applicable to bare n-particle interaction 
vertices. The self-energy is equal to the one-particle vertex, 
£= Yn=l- 

A derivation of the flow equations for the vertex functions 
within the generating functional approach in Keldysh formal- 
ism is given in Ref. |38l Reference 25 describes an equiva- 
lent derivation of the flow equations based on diagrams and 
provides diagrammatic rules for their formulation. The result- 
ing flow equations form an infinite coupled hierarchy where 
the flow of the n-particle vertex functions, dy^/dA, is a func- 
tional of , 72 , . . . , 7^ + i- For example the flow equations for 
the one- and two-particle vertex function read 



dA 



'l' 1 



1 

2k 



and 



dA^'2'|12- 



2^^2'3'|123 ^3|3' 
+ 2K ^"'2'|34 5 3|3' G 4|4' T^'4'\12 



^rf'4'|32 



2^3'|14 



c^X s~~<X 1 s-<X oX 
^313' 414' ~ t ~ 313' ^414' 



c^X s~iX I s-<X oX 

^313' 414' ~ t ~ 313' ^414' 



(27) 

(28a) 
(28b) 

)^2'|14 

(28c) 

$2'|32- 

(28d) 



and indices occurring twice in a product implicate summation 
over state and Keldysh indices and integration over indepen- 
dent frequencies. Furthermore, 



G x g x l s x g x x G x 



with 



dA 



(30) 



(31) 



denotes the so called single scale propagator. We call the 
contributions d28bb . d28cl) , and d28db to the flow of yf/ 2 /|i 2 
particle-particle, exchange particle-hole, and direct particle- 
hole channel, respectively. 

For practical computations the exact infinite set of flow 
equations is reduced to a closed finite set of approximated flow 
equations. Typical approximations are to neglect the flow of 
higher order vertex functions by setting dy^ /dA = for n > «o 
and to ascribe an effective parametrization to the remaining 
ones. In the present investigation of the SIAM we neglect the 
flow of . Two different parametrizations of ^ will be dis- 
cussed, a static (i.e. frequency independent) and a dynamic 
(frequency dependent) one. 



IV. FLOW PARAMETER 

In Ref. |42| it has been shown how the choice of the flow 
parameter determines whether certain exact properties of the 
vertex functions are conserved by an approximated fRG flow. 
For the study of the SIAM the conservation not only of causal- 
ity but also of the Kubo-Martin-Schwinger (KMS) conditions 
characterizing thermal equilibrium is important. The correct 
description of thermal equilibrium is necessary in order to 
capture fundamental aspects of the low energy properties of 
the model. Consider for example the expansion of the self- 
energy at eVg = ji = 0, B = in leading order in frequency, 
temperature, and voltage, 



E^(c»,r,V)~|+(l-£±£| ( o 



(2s Xc) 

1 4T 



a 2 + {KTf + heVf 



(32) 



where Xc and Xs denote the reduced charge and spin sus- 
ceptibility 47 - 62 . The Fermi liquid behavior[32l determines the 
proper shape of the spectral function, in particular its height 
and width. The derivation of ( 1321 ) is based on thermal equilib- 
rium particle statistics and hence requires the validity of the 
KMS conditions for a treatment within Keldysh formalism. 

Hybridization can be used as a flow parameter which con- 
serves causality and the KMS conditions 42 . For that purpose 
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the constants Tl.r and T from Eqs. (fl~8T > and ( fT9l are enhanced 
artificially by setting 



T - 1 

r ( /> = r,. + -^A, r = L,R, 



(33a) 
(33b) 



where A flows from °° to 0. Via the hybridization, the com- 
ponents of the noninteracting reservoir dressed propagator ac- 
quire the A -dependence 



fe) " t(co) = tH - £g+ i 1 (r+A)/2 ' 



(34a) 



fcl)£ V («) = tel)a e W> (34b) 

(jg X )a(a»)= [l-2/ eff (0))] [(rf>)-feO£ V («>)], 

(34c) 

where /^(o) = / e ff(d)) does not depend on A since 

rl^/r^ = r r /r, compare Eq. (fJTJ. The components of the 
free single scale propagator^ = dgi/dX are given by 

sV(v) = - i -[ 8 r(co)]\ (35a) 

4 v (ffl)=^ et (fi))t (35b) 

£(a>) = [l -2/ eff (<»)] [#» -*f (»)] . (35c) 

and those of the full single scale propagator Sx = 
Gxg^ng^Gx by 



Sf(fl») =5f(o) t , 

(o>) = - l -Gf \co)Gf(co) + ^Gf(cD)G{™(co) 



(36a) 
(36b) 



-i[l -2/ eff («)]Gf(0))Gf (to). (36c) 

A special situation occurs when 

G» = [l-2/ eff (a))] [Gf(0))-Gf (*>)]. (37) 

Then Eq. ( 136cl i can be simplified to 

Sf(co) = [1 -2/ eff (o>)] [sf '(to) -Sf (fl>)] . (38) 

The condition ( 1371 ) is fulfilled for instance, due to the fluc- 
tuation dissipation theorem, in thermal equilibrium, when 
jU L = jU R and f L ((0) = /r((o) = f{a>) = f e «(co). In non- 
equilibrium inelastic interaction processes mediated by the 
two-particle interaction will in general break the relation (f37l >: 
via contributions to E K and to the anti-Hermitian part of E Ret 
they tend to smoothen the effective distribution. Later we will 
feed back only the static part of the self-energy into the RG 
flow so that (|37T i is fulfilled automatically, see section lVI.DI 
The starting values of the vertex functions at A — °° can 

be determined as follows. For A — > °° the propagator g?' a 

vanishes as 1/A whereas fdcog^" (<X>) approaches a finite 
constant. Hence all diagrams having more internal lines than 



integrations over independent frequencies vanish. For a di- 
agram with m two-particle vertices which contributes to the 
n-particle vertex function y„ the number of internal lines is 
2m — n while the number of integrations over independent fre- 
quencies is m — n + 1. Therefore only the first order contribu- 
tions to the vertex functions do not vanish for A — > °°. This can 
also be understood in terms of time dependent diagrammatics: 
Due to the diverging decay rate X —> °° only diagrams which 
are completely local in time do not vanish. These are exactly 
the first order diagrams, namely the Hartree-Fock diagram for 
the self-energy E and the bare interaction vertex for y 2 , 



v Ret,Av / 
h a (A 



> = o, 

)=v, 

) = lim 

_ U 
~ 2 

E K (A =00) =0. 
In Eq. d39b| ), v is given by 



n > 3, 



(39a) 
(39b) 



(39c) 
(39d) 



_a\ OL 1 \OL\ On 
C7 G\ <J2 



with 



2 v CT 1 't7^|(J 1 a2' 

0, 



if a[ + + a\ + a-i is odd, 
else, 



(40) 



a. a, I a\ 02 



U, if C7j' = G\ = & 2 = <J2, 

= {-U, ifa[ -ai = a 2 -cj 2 , (41) 
else, 



compare d25l l and (HJ. The integral in (I39ct has been evaluated 
by 



1 f ±x M£q'+xT x /2) 



% J- 



dx- 



1 



1 

2' 



(42) 



where x = 2 ( (0 — e a i ) /T^ . The physical interpretation of d42l) 
is that the mean occupation of an infinitely broadened level is 
1/2. 

Hybridization as flow parameter has the advantage to con- 
serve causality and the KMS conditions (the latter being vio- 
lated for example by the imaginary frequency cut-off intro- 
duced in 25) and to be applicable to zero-dimensional systems 
(as opposed to the momentum cut-off). Since it changes the 
free propagator in a smooth way, Sx is not restricted sharply 
to a single scale by a delta function, as it happens for step- 
function cut-offs, e.g. the momentum cut-off. Such a delta 
function simplifies the flow equations since it cancels one of 
the necessary integrations. The hybridization-flow thus has 
the disadvantage of an extra integration in the flow equations 
compared to step function cut-offs. 
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V. FRG IN STATIC APPROXIMATION 

Before discussing the more elaborate frequency dependent 
truncation scheme, we describe a basic static approximation to 
the flow: the flow of the n-particle vertex functions for n > 3 
is neglected and the flow of the two-particle vertex function is 
reduced to the flow of a frequency-independent effective in- 
teraction strength. As a consequence the self-energy remains 
frequency independent which means that the influence of the 
interaction on single-particle properties of the system is de- 
scribed by a mere shift of the single-particle levels. This ap- 
proximation has been studied at T = within a fRG treatment 
based on Matsubara formalism using an imaginary frequency 
cut-off 54 - 55 . There it was shown to produce a Kondo scale ex- 
ponentially small in U /T which occurs in the pinning of the 
renormalized level position to the chemical potential. Results 
for the dependence of the linear conductance on the gate volt- 
age for different magnetic fields have been found to be in very 
good agreement with Bethe-Ansatz and numerical renormal- 
ization group computations. It turns out that in thermal equi- 
librium at T = the hybridization flow produces the identical 
flow equations. Hence this Keldysh approach incorporates all 
features found in Ref. 55 for the Matsubara fRG. 

Apart from its success in describing the linear conductance 
the applicability of the approximation in question is restricted. 
Given a frequency independent self-energy, the spectral den- 
sity Per (ft)) is a Lorentzian that is centered at the renormalized 
level position e a and has fixed width and height, 



1 



r/2 



n((o-e a ) 2 + F 2 /4 



(43) 



Therefore the approximation cannot describe any features 
connected to details of the spectral function such as the forma- 
tion of a Kondo resonance with side bands and its dependence 
on temperature or voltage. That's why we do not elaborate 
in detail on this approximation but merely show that at zero 
temperature and in equilibrium the resulting flow equation are 
identical to those of Ref.l55l 

Reducing the two-particle vertex function to a frequency 
independent effective interaction means setting 



/_ x cc j a 2 1 «i «2 



(44) 



where \>x is obtained from v by replacing the bare interaction 
U by a renormalized one Ux in ( HTt . Here the effective inter- 
action Ux is a real number whose flow starts at 



U x=m . = U. 



(45) 



real, so that we can speak of a flowing effective level position 
= e a ^ ne initial condition for the level position is 

according to Eqs. <£3j and ( I39cb 



E**(A=oo)= e y g -fffl. 



(46) 



The starting values (l45T l and d46l) are identical to those of 
Ref. |55j Note that in the framework used there the initial con- 
dition for the level position is the result of a first stage of flow 
from A = °o to A = Ao — > °°, see e.g. Ref.l54l 

From Eq. (ffTJ) we derive the flow equation for the level po- 
sition 



de, 



_£L = 

dA 2n 2 

Ux 



dcoSf w (co) 
g/dfi)sign(fi)){[G^(o))] : 



where we used Eq. (138b with 

1 - 2/eff(G>) = 1 - 2/(0)) = sign(o) for 
Inserting 



[G^(O))] 2 }, 
(47) 



0, n = 0. 
(48) 



1 



Gf(o)) t 



fi)-e*+i(r + A)/2 
we can evaluate the integral in Eq. ( l47l ) and obtain 



(49) 



de£ Ux 



£- 



dA 2k £ i 2 + ( y + x y/4 



(50) 



The flow equation (1281 1 for ji generates a frequency depen- 
dent vertex function which has a richer structure in terms of 
Keldysh and spin indices than indicated in Eq. d44l >. It can be 
shown that after neglecting the contribution from the three- 
particle vertex function (12 8 at this structure can be reduced 
to the simple form given in Eq. d44b by setting all outer fre- 
quency arguments equal to the chemical potential fi. We note 
that this projection onto the Fermi surface is the core idea of 
the static approximation. The resulting flow equation acquires 
the form 



dA 



(/f)^ 12 (0) + (/f)^ |21 (0) 



+(/f)?l |22 (o) + (TO"(o) 



ph\22| 12 / 



(51) 



As a consequence of that approximation to the vertex function, 
the flowing self-energy remains frequency independent and 



where we abbreviated frequency integrals appearing in the 
particle-particle and particle-hole channel by 
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(/r)S aK (») = ^/d»' 



a i i «i 



(52a) 
(52b) 



Making use of Eqs. OfJal) . ( f36bl , (O, OH), 6D and @9) we 
evaluate the integrals and obtain 



dA 



U 



(r + A) 2 /4 ef +(r + A) 2 /4 



(53) 



The substitution A — > 2A maps the flow equations ( f50b 
and d53l onto those of Ref. HH We note in passing, that the 
static approximation scheme can be carried out as well with 
the imaginary frequency cut-off from Ref. [25| instead of the 
hybridization flow parameter; the special advantage of the 
hybridization flow to conserve the KMS relations becomes 
relevant only in the dynamic approximation scheme intro- 
duced in section [VI] The flow equations of the static scheme 
based on the imaginary frequency cut-off are again identical 
to Eqs. d50l > and ( 1531 and to the flow equations of Ref.HH The 
real frequency cut-off used in Ref. 38 on the other hand pro- 
duces different flow equations which nevertheless lead again 
to the same final solution^. 

It is easy to generalize the static approximation scheme to 
non-equilibrium. Since the Keldysh component of the self- 
energy is not renormalized we merely need to replace d48| ) by 

1 - 2/ eff (£») = £ — sign(o) - jU r ) , (54) 

r 

which leads to a superposition of the flow equations found for 
thermal equilibrium in the form 



d4 = u ly r L 

dA 2n^ r (4-aO 2 + (r + A) 2 /4' 

dA t r ,.l .. ^2 



(55a) 



* r r ( e A_ Mr )2 + (r+A)2/4 



(55b) 



This approach, however, does not provide a reliable descrip- 
tion of the influence of bias voltage, since it restricts the spec- 
tral function by construction to a single Lorentzian peak; a 
splitting or damping of the peak is principally impossible. 



VI. FRG IN DYNAMIC APPROXIMATION 

In order to overcome the restrictions of the approximation 
described in section [V] a frequency dependent self-energy is 
required. The flow equation d27T i produces a frequency de- 
pendent self-energy only if the two-particle vertex function 



is frequency dependent. Therefore our aim is to extend the 
scheme described in section [V] in a way that a frequency de- 
pendent two-particle vertex function is generated. We seek 
for an extension which is as basic as possible; in particular we 
neglect the influence of the three-particle vertex function fur- 
ther on, such that the flow of the two-particle vertex function 
is induced by the three contributions ( |28b4 - |28cfl >. 



A. Ladder approximations 

Recent investigations of the SIAM by diagrammatic tech- 
niques have revealed the importance of the exchange particle- 
hole ladder (RPA series) for the emergence of Kondo 
physic s 75 i 76 . Let us therefore in a first step neglect all contri- 
butions to the flow of the two-particle vertex function except 
for the exchange particle-hole channel (|28ct . If we suppress 
additionally the feedback of the self energy into the propa- 
gators except for the initial Hartree term J39cb then the flow 
equation for the two-particle vertex functions reads 



dA 



(7x) 



1'2'|12 



2% 



(tx) 



1'4'|32 



x (Yx. 



3'2'|14 



'33' 



dA 



-(fx) 



3'2'|14i 



(56) 



where is obtained from Eq. d34i l by replacing e a with 
£ CT + U /2. We denoted the vertex function in this approxima- 
tion with a superscript "x" which refers to "exchange". The 
solution of the flow equation (l56l > with the initial condition 
y£ =00 = v is given by the exchange particle-hole ladder (RPA 
series) 

(7a)i'2'|12 = Vl'2'|12 + ^" f7 l'4'|32^3|3'^4|4'' f '3'2'|14 + ■ ■ ■ 

= 'f 7 l'2'|12 + ^'l'l'4'|32^3|3'^4|4'(7A X )3'2'|14- (57) 

Due to frequency conservation at each interaction vertex, the 
frequency dependence of the ladder is reduced to a single 
bosonic combination of the external frequencies, 



Y Z( (a [M\au(*2) = (YD(X), 



with 



X = a>2 - a>i = a>2 - <*>[ ■ 
We evaluate the RPA series and obtain 



(58) 
(59) 



,.,12122 , V ,_U ,U 2 (g|WPQ ,, m 

05 W^ - y + Ti-f (fi j) ff ,(x)' (60) 
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where 



equations satisfy 



+ *?a(® "§)«&(<» + §)] (6D 

denotes the particle-hole polarization operator. In the special 
case T = 0, B = 0, V = 0, eV g = /X = we find 



nX(X-T x ) 



•In 1 



r772 



(62) 



Expanding 1 / (B\) a -^{X) in powers of X/T^ results in 



2" 



CP 



(63) 



27rX-i(r A /2-f//7r) 
which features a singularity when A reaches the value 

A = A c =Wj%-T. (64) 

12122 

Additionally, for A < A c the function (y^) a l^ a —(X) exhibits 
the wrong analytic behaviour, being analytic in the upper half 
plane of X instead of the lower one as required by causality 42 . 
Therefore the flow A — > can only be finished if 



U <U C = izT/2. 



(65) 



This limitation is not even overcome when the full self energy 
is fed back into the flow, replacing s and g in Eq. (l56*i l by S 
and G. 

Hence we also have to take into account the other two chan- 
nels d28bl28dt . A parquet summation based procedure mix- 
ing the exchange particle-hole and the particle-particle chan- 
nel has been set up in Ref. 75, in which the authors study the 
SIAM within Matsubara formalism. Essentially the particle- 
particle channel serves them to renormalize the interaction U 
which enters the exchange channel to an effective value U e s 
which is always lesser than U c . The difference U c — U e g is 
then used as a measure for the Kondo scale. 

We aim to implement a similar proceeding within the 
framework of the functional RG. We consider it appropriate to 
include also contributions from the direct particle-hole chan- 
nel, because this channel bears the same singularity as the ex- 
change channel. In order to see this let us define the particle- 
particle and the direct particle-hole ladder on the analogy of 
Eq. ( |56l l by the flow equations 



d 

dl 



TT-(7a)i'2'|12 



i , p. d 4|3 



dA 

A 

313' 



A J3'4'|12) 



<7i 



(66a) 



(66b) 



dA ^ 4 ' 2 'l 32 ' 
In correspondence with Eq. d57l ) the solutions of these flow 



(7;£)l'2'|12 — Vl'2'| 12 + ^—Vi'2'|34g3|3'£4|4' ( 7^ ) 3'4' 1 1 2 , 
(7a)i'2'|12 = Vic 2 '| 12 — IT- Vl'3'|14^3|3'^4|4' (7a)4'2'|32- 



Their frequency dependence takes the form 

{yl){co' u c^\co u ay 1 ) = {yl){n) 
{y*){co[,ai\co u a> 1 ) = {y*){A) 

with 



II = coi + a>2 = <x>[ + (0^, 
A = co[ - coi = (O2 - a>2- 

Evaluating (|67a| ) one finds 

. PU2|22 (m _U , t/ 2 (^)tfff(Il) 



2 4 l-|(BP)^(n)' 



(67a) 
(67b) 

(68a) 
(68b) 



(69a) 
(69b) 



(70) 



where at T = 0, B = 0, V = 0, eV g = = the particle-particle 
polarization operator B p ^ takes the form 



1 ^ infi-i n 

a:n(n+ir A ) V ^1/2 



(71) 



such that 7^ is a regular function on the real II-axis for all 
values oft/. In contrast, (7^) 12 ' 22 satisfies 



2 

1 r 



(72) 
(73) 



and thus possesses the same singularity as the exchange 
particle-hole channel. 



B. Approximated mixing of the channels 

In an approximation taking into account all three channels 
(I28bl428d"l ). the two-particle vertex function will depend on all 
three frequencies, y(U,X,A). Due to frequency conservation 
these three ones are indeed sufficient to express the general 
frequency dependence. In the following the three bosonic 
frequency arguments of the vertex function will always be 
indicated in the order (II, X, A). When all three contribu- 
tions ( T28bU28dl i to the flow are taken into account simulta- 
neously, then the dependence on II, X, A is mixed through the 
feedback of the vertex function on the right hand side of the 
flow equation that reads 
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d 

dX 



ri 2 , ]12 (n,x, a) = ±-J do)|^ 2 ,| 34 (n, (»+ 2^, © - 2^) s* |y (§ - ©) g^, (n + to) $ 4 , |12 (n, + a,, i+a _ 

(74a) 

- ^132 + S±A - 0,) [S* l3 , (ft) - |) G£ |4 , (CO + f ) + G*, 3 , (tO - f ) S 4 A |4 , (tO + f ) ] ^ 14 (fl> + 0) - 2=A) 

(74b) 



A'2' 32 I 2 + 



to 



n+x 

' 2 



to, A) 
(74c) 



r 



Here we use shorthand notation like 

7rai2= (tf) ' 



a j | a { a 2 



(75) 



and indices occurring twice in a product implicate summa- 
tion over state and Keldysh indices. The flow equation (l74l 
leads to a complicated dependence of the vertex function on 
the three frequencies (II, X, A). As a consequence a numeri- 
cal solution of the flow equation requires a sampling of three 
dimensional frequency space which constitutes a high com- 
putational effort. The structure of the flow equation suggests 
to approximate the frequency dependence of the two-particle 
vertex function by 



7 A(n,X,A)~v+(pP(n) 



-<pl(X) + (pi(A), (76) 



where the bare interaction vertex v is the initial value at the 
beginning of the flow, v = Yx=°°, and where (p p (YY), qjf (X), 
9? (A) are approximations to the parts of y produced by the 
three channels J74ab . d74bt . (T74cl >. respectively. An analo- 
gous approximation has been investigated in studies of the 
SIAM based on the equilibrium Matsubara fRG 57 . The re- 
sults achieved with and without this approximation were of 
the same quality. Reliable results were obtained for small and 
intermediate U /Y. The limiting effect for large U /T seems 
to stem from omitting the three-particle vertex function. For 
small U /Y the system can be described by second order per- 
turbation theory which also produces a two-particle vertex 
function of the form d76i l. Our fRG treatment will fully com- 
prise the diagrams of second order perturbation theory for the 
self-energy and for the vertex function. Therefore, for U — > 
the fRG results asymptotically approach those of second order 
perturbation theory. 

In order to achieve the form (l76l l where the frequency de- 
pendence is split up into three functions we have to elimi- 
nate X and A from d74aK n and A from (|74b| i, and II and 
X from ( 174-cb . It is obvious that any manipulation of this 
type can only yield convincing results if the frequency depen- 
dence of the vertex functions is not very pronounced. This 
limits the range of applicability of the approximation to small 
and intermediate interaction strengths. Kondo physics emerg- 
ing for large U /Y cannot be described in general: we expect 
Y(Yl,X, A) to exhibit a sharp resonance as function of X in this 
regime (cf. Ref. 75 and section [VI. At ; furthermore omitting 
the three particle vertex is not justified for large U /Y. 



Different approaches in order to get rid of the frequency 
mixing on the right hand side of d74i i are conceivable. The va- 
riety of reasonable replacements is however restricted by the 
necessity that the components of the vertex function maintain 
their fundamental features concerning exchange of particles, 
causality and the KMS conditions as described in Ref.H^- We 
have tested different possibilities. Here we present only that 
one which provided the most convincing results. It is based 
on a procedure to assign a static vertex to the functions <pP' x ' , 



9l(n) 

9f(A)->*i- 



(77a) 
(77b) 
(77c) 



The < I ) P' x ' d are frequency independent and represent effec- 
tive static interactions between particles of either opposite or 
identical spin state. Since permutations of the incoming in- 
dices map the exchange particle-hole channel onto the direct 
particle-hole channel and vice versa, the <J> x d have to satisfy 



(<&; 



x , a[ a 2 1 a 2 a.{ 
X / ctct | a a 

x , a[ a 2 1 0C[ a 2 

X 'aa\aa 
x , aj a' 2 1 0C[ a 2 

X I aa\aa 



^ X / (TO 7 ! cf CT ' 



(j ,a[a 2 \a 2 ai 
X i cj(7 1 aa '■ 



\ X'aalaa 



(78a) 
(78b) 
(78c) 



This is achieved in the parameterization 



I <J\ (7? 





if a[ = 


0\ = 


°'2 = 


<?2, 


-u p 
u x> 


if a[ = 


0\ = 


°2 = 


<*2, 


0, 


else, 








u x 

U V 


if a[ 


= 0\ 


= a' 2 


= ^2, 


-f/ d 


ifa[ 


= o x 




= C*2, 




, if a[ 


= CJ| 


= o> 2 


= o 2 , 


0, 


else, 










if a[ = 


0\ = 


o' 2 = 




—TI X 
U X> 


if a[ = 


o\ = 


o' 2 = 




VV X CT] ' 


if a[ = 




G> = 


<?2, 


0, 


else, 









(79a) 



(79b) 



(79c) 
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and 




if a[ 
else, 



a\ + (X2 is odd 



(80) 



with t/P' x,d and Wf being real numbers, compare (l40l 
and ( 1411 ). The detailed procedure to determine appropriate 
constants , W d CT has to be consistent with the flow equa- 
tions and is discussed below in section lVI.CI In J79a[ > we ex- 
clude the possibility that an effective interaction between two 
particles of identical spin state could result from the particle- 
particle channel. The reason is that a static interaction vertex 
between particles of identical spin state vanishes when being 
anti-symmetrized since the direct and exchange terms cancel 
each other. This is consistent with the fact that such a con- 
tribution is not generated by the flow equation ( 181 at below, 
since v + + <f> d allows only interactions of particles in op- 
posite spin states. Also the fact that E/?' x,d are independent of 
of the spin (while W d is not) is consistent with the structure 
of an effective interaction vertex, compare (14 It . and will be 
reproduced by the flow equations below. 

When inserting now the form d76l ) on the right hand side of 
the flow equation d74l ), the sum structure is reproduced by the 
flow if we replace cpf(X) and <p d (A) in the particle -particle 
channel ( |74at by the constant vertices 4>^, <t> d and proceed 
analogously for the other channels: in the flow of any channel 
the other two channels are reduced to constant vertices. The 
flow equations for <p^' x ' d are then given by 



^K)i' 2 '|i2(n) = i(v+(pP(n)+4> x ,+4> d ) 



1'2'|34 



(if ) 34 |3'4< (n) (v + cpl (u) + <t>\ + 4>i) 3 , 4 , |12 , 

(81a) 
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x(lf) }4l y 4 ,(X)(v + <Pl + Cpl(X)+^) y2 , ll4 , 

(81b) 

^(<PA d )i'2'|i2(A) = -(v + 4>P+*I + (p d (A)) 1 ,3, |14 

x (/P h ) 34 | 3 , 4 ,(A)(v + 4>P + <t>\ + (p d (A)) 4 , 2;|32 , 

(81c) 

where we used a summation convention for index numbers 
representing the state a and the Keldysh index a, and 



/rPP.ph\ , , , 
V l X ]34|3'4' 



/■rPP.ph-v«3«tl«3«4 

VX J0304 



(82) 



refers to d52l l. 

The different spin and Keldysh components of the functions 
(pPAd are no t completely independent of each other: there ex- 
ist identities connecting them which stem from the invariance 



under exchange of particles and complex conjugation . Fur- 
ther identities follow from spin conservation and from the spe- 
cial structure of the flow equations (ISTI ). In Appendix lAl we 



determine a set of independent spin and Keldysh components 
of the functions (pP x d from which all other components can 
be derived. It is then sufficient to compute the RG flow of 
these independent components; the corresponding flow equa- 
tions are presented in the Appendices iBl and ICl 



How to determine Uf x ' A and 

A a a 



The approximation scheme described above requires a pro- 
cedure which assigns effective constant interaction vertices 
4>P ,x ' d to the functions <p£(n), cp\(X), cp^(A). Hence it is 
a central question how to determine appropriate real num- 
bers U? , Wjj 1 which characterize the vertices <$^ x ' in 
Eq. d79l ). We propose a very simple scheme for the case of 
thermal equilibrium. In order to generalize this approach to 
non-equilibrium however we will need to formulate additional 
constraints. 

Let us first assume thermal equilibrium. Then we can show 
that 



*5 

n 



<Pl(X = 0), 
<J»d =( pd(A = 0) 



(83a) 
(83b) 
(83c) 



have exactly the structure given in Eqs. (T79b and dSOt . For the 
proof we make use of the special spin and Keldysh structure 
of (pP' x ' d that is described in Appendix lAl 

We first discuss the particle-hole channel. From 
Eqs. dA15ab and (lA18ab we infer that (a d ) CT tf(0) and 
(fl d )(jct (0) are real numbers. Using the KMS conditions (IA16I ) 
and (|A19b combined with Eq. dC8l ) we conclude further that 
(^ d )(7o : (0) = and (b^) aa (0) = 0. In thermal equilibrium we 
can use the fluctuation dissipation theorem 

Gf(co) = [1 - 2/(0))] [Gf\(D) - Gf (co)] , (84a) 

Sf(a>) = [1-2/(0))] -Sf (co)} , (84b) 

to show that 

(^ ph ) 2 i 22 (o) + (/f) 2 i 12 (o) = 

-2Re|^|do) [G£(co)Sg(co)+S^(co)G£(co)] 

x[l-2/(o))]l (85) 



is a real number. From the flow equation (IC6ab it follows that 
(a\) a s(0) is real. The KMS condition JA13t then entails that 
(b\) a o(0) = 0. In total this means that we can set 



4> d =(p d (0) and ^=<Pa(0), 



with 



ut 



w 



X a 



2(4)aa(0), 

= 2(^)^(0), 
2(ai) aa (0). 



(86) 

(87a) 
(87b) 
(87c) 
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When we insert ( I87cb into (IC7ab this yields in particular 
d(a^) C7 o , /dA = 0. Combining this with the initial condition 
(a^ =0O )<TCT = we conclude from d87bb that 



W = 0. 



(88) 



For the particle-particle channel we combine the fluctuation 
dissipation theorem d84t with 



l-2/(fl)) = -[l-2/(-fl))] 
(mind the convention /I = 0) to find that 

(/f)S 1 >) + (4 pp )S 21 (o) = 



(89) 



2Re{^/dc 



cr(«)C(-«)+5r(«)G# v (-«) 



Ret/ 



-.Av i 



x [1-2/(0))] | 



(90) 



is a real number. Hence it follows from the flow equa- 
tion ( lC3al > that (a P x ) a o{0) is real. The KMS condition (lAlOb 
then demands that (#! = 0. Therefore we can set 



with 



^1 



; 2(4)^(0). 



(91) 



(92) 



(Note that ( v + + <i>\ + 4> d ) = is exactly the flowing 
effective interaction vertex used in the static approximation 
scheme in section [VI ) 

In the case of non-equilibrium with fi^ + /Ir = we will 
also use d87} and d92}. While Uf and W A d are then again guar- 
anteed to be real valued by ( |A15at and ( |A18at , additional 
assumptions are necessary to ensure this for U£ and Iff. If 
the propagator which we feed back into the flow of the two- 
particle vertex satisfies Eq. d37| i, then Eq. ( f85T > is still valid in 
non-equilibrium, with / e ff(ft)) replacing /(ft)); this warrants 
that I/£ is real. Under the additional assumptions 



T L = T R and r L = J R 



(93) 



also Eq. d90t can be maintained; in that case also IT? is real. 



While this identity is guaranteed in equilibrium by the 
KMS conditions, it constitutes an approximation in non- 
equilibrium; interaction effects that lead to a redistribution of 
particle statistics are neglected. 

The results discussed in section IVIIII below show that the 
feedback of the full frequency dependent self-energy into the 
flow leads to an artificial smoothening of spectral features. 
The Kondo resonance at 0) = for instance acquires a much 
broader shape than expected. Better results are achieved with 
a different feedback scheme for the self-energy which reduces 
its back-coupling into the flow to a simple shift of the single- 
particle levels. In analogy to the way of determining an ef- 
fective static interaction vertex from the vertex function in 
section IVI.CI we assign a real and static level shift E to the 
self-energy by setting 



E* ea = E£ VA = ReE 



Av A 



-Ret X 



(ft> = 0), 



7 Kl 



0. 



(95a) 
(95b) 



Only this static renormalization enters the propagators on the 
right hand side of the flow equations for T.^ and . Due to 
the missing Keldysh component of the condition ( f37b is 
then readily fulfilled. 

It turns out that this static self-energy feedback produces re- 
sults distinctly better than in case of the full frequency depen- 
dent self-energy feedback. We also studied mixed schemes, 
for instance, to feed the static self-energy into the flow of 
while the full self-energy is used for the flow of E^, or vice 
versa. However, static feedback into both flow equations per- 
formed better. Hence all results shown in section IVIIII are 
computed with this scheme. The fact that static feedback per- 
forms best highlights the difference between the Keldysh fRG 
and its Matsubara counterpart, where the full self-energy feed- 
back is used 57 . Note that in frequency-dependent approxima- 
tions the two methods are not equivalent (in equilibrium), even 
if the same feedback scheme is used. We suspect that the feed- 
back of the full self-energy into the flow becomes favorable in 
a truncation scheme taking into account contributions from 
the flow of the three-particle vertex. Such a scheme, allowing 
for a renormalization of the static part of the two-particle ver- 
tex, might lead to an enhanced effective interaction strength, 
higher values of the effective mass, and a sharpening of spec- 
tral features - following the tendency expected from the Ward 
identities 4 -"^. 



D. Self -energy feedback 



VII. APPROXIMATE SOLUTION IN THE 
PARTICLE-HOLE SYMMETRIC POINT 



The propagators S, G appearing in the flow equation for the 
self-energy and for the vertex function are full propagators, 
which means that they are dressed with the current value of the 
self-energy E^ (ft)). In this way the self-energy feeds back into 
its own flow and into the flow of the vertex function. In the 
preceeding section we described that our flow scheme requires 
the validity of the fluctuation dissipation theorem Eq. d37l i. or, 
equivalently, 



In this section we discuss the special case of vanishing mag- 
netic field, B = 0, and particle-hole symmetry; the latter is 
given for eV g — (/II + A*r)/2 = 0. Note that we already in- 
troduced before the restriction Tl = Tr. In that symmetric 
situation the RG flow does not generate any contribution to 
ReE^ et (o) = 0) which hence remains at its initial value U /2 
given in Eq. (I39cb . The renormalized effective level position 
is then given by 



E£(ft)) = [l-2/ eff (ft))] [Ef»-Ef(ft)) 



(94) 



e^=e CT +ReEfJ(ft) = 0) = 



(96) 
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compare Eqs. (01 and (195 at . Hence the effective single- 
particle levels for the two spin-states are degenerate and fixed 
in the middle of the two chemical potentials throughout the 
flow. 

This section is organized as follows. At first, we discuss 
the structure of the vertex functions in the particle-particle 
(PP) and exchange particle-hole (PH(e)) channels. At sec- 
ond, we focus on the vertex functions in the direct particle- 
hole (PH(d)) channel. We distinguish this case from the other 
two channels because of a bit more involved spin structure. 
At third, we formulate a simplified version of the flow equa- 
tion (f27T > for the self-energy. On its basis we recover the 
Fermi-liquid relation s 47 ' 62 for the (imaginary part of) self- 
energy providing an analytical fRG estimate for a character- 
istic energy scale of the quasiparticle resonance at zero fre- 
quency. 



A. Particle-particle and exchange particle-hole channels 

Let us first discuss the structure of the vertex functions in 
the particle-particle and exchange particle-hole channels char- 
acterized by energy exchange frequencies II and X, respec- 
tively. 

From Eqs. d87bb and d88l we establish the property 



where the static part 



U x = U + U{ + Ul = 2a x (0) = 2a 3 (0) 



obeys the equation 



(108) 



(109) 



In case of particle-hole symmetry and zero magnetic field the 
following property holds 



F 3 (X) = -F l (-X) 1 



(HO) 



which, in particular, implies F\ (0) + Ft, (0) = 0. Therefore the 
static component Ox = U does not flow within the approxi- 
mation introduced in subsection [VLB] 

The flow equations ( IC3bb and ( IC6bb for the functions j 100b 
and ( 11041 > read in the new notations 



dfci(n) 

dA 

db 3 {X) 
dX 



= | fll (n)| 2 //!(n) 

+ 2iim{fli(n)6i(n)Fi(n)}, (in) 

= \a 3 (X)\ 2 H 3 {X) 

+ 2ilm{a 3 (X)b 3 (X)F 3 (X)} . (112) 



(«iW(0)=^/2 = 0, 



(97) 



which leads, as we will show below, to decoupling of PP and 
PH(e) channels from PH(d) channel. 

Taking into account d97l i we introduce for convenience the 
following shorthand notations (cf. Eqs. ( f79l , JA9t and ( lA12b ) 



fli(n) 

Fx (II) 

fti(n) 
#i(n) 



(4W(n) + ^i 

(ID 



'X 
2ilm 



f)S n (n) + 5(/f)S 2 ( 



and 



a 3 (X) = 

F 3 (X) = 

b 3 {X) = 

H 3 {X) = 



u +u? 

—A + (al) a c(X), 



2ilm 



(98) 

(99) 
(100) 

(101) 



(102) 

(103) 
(104) 

,(105) 



where the functions ( I98l l and (I102l i are specifically defined to 
fulfill fli(0) = 03(0). The flow equations ( lC3ab and (lC6ab can 
be then transformed into 



dai(n) 
dA 

da 3 (X) 



ai(n)] 2 Fi(n) 

Ox- : 



F 3 (0), (106) 



dA = It' F M+w)] 2 m), do?) 



The functions F\ t3 and Hi 3 are given by 

F(n) 



-G^ + a>>)G^-co>) 



(113) 



F 3 (X) 



GT(co'~\)GUco' + \) 



+G$(cD'-|)G^(fl)' + |) 



(114) 



and 



ffj(n) 



dco' 



J__d_ 
2n~dX 



(115) 



H ^ X) = Tnlk! d(0 ' 

-.Ret ^Av 



x [GTGp + G^G^ + G^G§\, (116) 

where the partial derivative d/dX is supposed to act only on 
explicitly A -dependent term within G's, and the arguments 
of Green functions in ( II 15l l and ( II 16b correspond to those 
in ( II 13b and ( II 14l i. respectively. (For brevity we also denote 
Gx = G.) By a straightforward evaluation one can check that 
in equilibrium the following identities hold 

ff!(n) = cofc^J[Fi(n)-z?(n)], (117) 

H 3 (X) = -coth(^\[F 3 (X)-F 3 *(X)}. (118) 
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With their help we can establish that the combinations 

Cl (n) = 6i(n)-cofh^J[ai(n)-flj(n)],(ii9) 



o(.Y! = />.,(* ) + colli ( ^ ) [^(.Y)-^(.Y)] 



obey the equations 



dci 

= 2Re{aiFi}ci, 
--7j- = 2Re{fl3F 3 }c3, 



(120) 

(121) 
(122) 



It explicitly manifests that the exponential scale develops in 
the broadening of the renormalized frequency-dependent ver- 
tex function in the exchange particle-hole channel. 



B. Direct particle-hole channel 

The flow equations for the vertex functions in this chan- 
nel are split off those considered in the previous subsection. 
We analogously introduce the new notations (cf. (IA15b and 
(EH) 



respectively. Since initially c\(Yl) = c 3 (X) = 0, they also re- 
main zero during the flow, accordingly to ( 11211 1 and d 1 22b . 
This proves that the KMS relations dAlOt and ( lA13t inherent 
to the equilibrium state are exactly preserved under the fRG 
flow in the chosen approximation. 

The derivative of the function F 3 at zero temperature reads 



dF 3 {X) 



dX 



x=o 



dco 



(123) 



x J 25(0)') [Gl et (co')lmG§{o)') -ImG^ et (o)')G# v (&)') 



020(A) 


= f + ( a l)aa{A), 


(129) 


a 2a (A) 


= ( fl i) ffff (A)--|£, 


(130) 




= (/f)S 21 (A) + (/P h )^ 22 (A), 


(131) 


b 2 o(A) 


= (bl)oo(A), 


(132) 


b 2 o (A) 


= (bi)aa(A), 


(133) 


H 2 „{A) 


= 2iIm|(/ ; [' h )^ 21 (A) + i(/P h ) 2 T 2 l 22 (A) 


K,134) 



+ h F ((o') 



Gl et (co')lm- 



jRet 

~lho> 



dG^(co') dG*«(a>>) 



-Im- 



dco' 



G^(co') 



and rewrite the corresponding Eqs. ( |C7t and dC9l ) 

d fl20 (A) 



where h F (co) = 1 — 2/(o)' 



sign (<b). Under constraints 



dA 



adopted in this section the formula ( 1123b simplifies to 



= -020(A) £ fl2 .v(A)F 2 ,(A) 

.v=U 



F 3 '(0) 



7T dX 



dco' 



+ 1^1 |Fi(0)+F 3 (0)], (135) 



45(0)') [lmGr(o)')] z + Ma)')^7 P^OO] 

(124) 



da 2(7 (A) 
dA 



-[a 20 (A)] 2 F m (A) - [a 2a (A)] 2 F 2a (A) 



^ ) ^(0), 



In the approximation achieved by neglecting the self-energy 
feedback Gx — > gx we can replace in d 1 14b the partial deriva- 
tive d /d A by the full one d/dA . Then we find that 



W) = dA 



%X{X-\Y X ) 



In 1+i 



X 



r A /2 



(125) 



is given by the full A -derivative of the particle-hole po- 
larization operator d62l . In particular, we obtain F 3 (Q) — 
—4/ [tzT 2 ^), and at small X < we approximately have 



db 20 (A) 
dA 



db 2a (A) 
dA 



(136) 



-a 20 (A)a*2c(A)H 2 a(A) - a* 20 (A)a 2a (A)H 2a (A) 

-a 20 (A)b m (A)F m (A) - a* 20 (A)b 2a (A)F 2 * a (A) 
-a* w (A)b 20 (A)F^(A) - a 2(T (A)Z> 20 (A)F 2cT (A), 

(137) 

-\a 20 (A)\ 2 H 2 c(A) - | fl2a (A)| 2 H 2a (A) 

-2Re{a 2a (A) J F 2(7 (A)}fo 2<J (A) 
-2iIm{^ (A)fc 20 (A)F^(A)}. (138) 



F 3 (X) 



^(2x-ir A ) 2 - 

This leads to an approximate solution of Eq. ( 11071 ) 

u _x_ ± mjM 

iK ' 2X-\B(U,Xy 

where the width of the Lorentzian is given by 

4U -M- 
B(U,X) = —e w i. 

K 



(126) 



(127) 



(128) 



We note that a 2 o(0) = 0^/2 and a 2cT (0) = as well as 
F 2a (A) — F^(—A) and H 2a (A) = H 3 (—A). One can also ob- 
serve that b 2 o(A) contains both real £> 20 (A) = Reb 2 o(A) and 
imaginary b l 20 (A) = Im/? 2 o(A) parts, while b 2a (A) is purely 
imaginary. 

Analogously to ( 111 8b the following identity holds in equi- 
librium 



H 2a (A) =coth 



—) 

2T J 



[F 2a (A)-F 2 : a (A)} 



(139) 
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It allows us to express the equations for b\ a (A) and for the 
combinations 

C2 o(A) = i4o(A) - coth ( [ fl20 (A) - o| (A)] (140) 



c 2CT (A) = fo 2CT (A)-coth^— ) [a 2CT (A)-aI CT (A)Kl41) 

in the following form 
db 



y 20 



dA 



dc 2 o 



= -R&{F 2a a 2a +F 2 aa 2 a}b 2 



20 



-iIm{F 2f jfl2(7 - F ST O2a}c20 
-iIm{F2 a a2o}c2(7 -i!m{F 2 aa2o}c2a, (142) 



= -iIm{F 2 CT«2<T - F2oaz5}b T 2o 



- Re{F 2a a 2a + ^2ct«2ct}c 2 o 
'Re{F2afl20}c2fJ -Re{F 2 afl20}c2a, (143) 



dc 2g 
dA 



2iIm{F 2 5a 20 + F 2<T a 2(T }fo2 



20 



~2Re{F 2 cra2o}c20 - 2Re{F 2a a 2a }c2a- (144) 



Since initially £> 2 o(^) = c 2o(A) = c 2o-(A) = 0, they remain 
zero during the flow [cf. Eq. ( 1A16H . This observation com- 
pletes the proof of the statement that the KMS relations are 
preserved in the chosen approximation scheme. 



C. Fermi-liquid relations for the self -energy 

In the case of particle-hole symmetry it is convenient to 
introduce the notations 



ai(n) = -ai(-n), 

d 2a (A) = -a 2cr (-A), 



(145) 
(146) 



and rewrite an equilibrium version of the Eq. ( IB lb for the re- 
tarded self-energy in a simplified form 



Ret 



dA E<J 



(co) 



dco 



U 



h F {co)lmsf l (co) 



+h F (co)di(cd - Co)lms^ -{-&>) 
+h B (co - co)lmd\(co - co)s^(-co) 
+h F (co)a 3 (cd - a>)Ims§- et (d)) 
-h B (cd - co)lma 3 (cd - co)s^ et (cd) 
+h F (co)d 2a (co- co)lms^\cd) 
-h B {co - co)lma 2a {cd - ct))s Ret (d))}j47) 

where h B (co) = coth ^t. In this equation the self-energy feed- 
back is neglected as well by replacing S% — > = Tjnf- 

Comparing Eqs. ( 1106b and ( 1107b with (1135b and (1136b . one 
can observe that the relation 



a 2a (co) 



a\(co) + a 3 (co) 



(148) 



holds at every A. The identities Ims Ret (co) = Ims Ret (— co), 
s Av (-C0) = -s Ret {co) are fulfilled as well. Therefore ( fT47b 
reduces to the form 



27Z 



dco 



x [h F (co)di(co - a>)Ims§- et (<») 

- h B (cd~co)lmdi(cd-co)s^ et (cd) 
+ h F (cd)a 3 (co- co)lms^ et (cd) 

- h B (cd-co)lma 3 {co-co)s^ et (cd)] (149) 
In particular, its imaginary part reads 

-^-ImE Ret (ct)) = ^- j dco\h F (cd) - h B (co - co)] (150) 
dA 2k J 

x lm[di{co- co)+a 3 (cd- co)]lms^ et (cd), 
which at T = amounts to 



AlmI-» 



(151) 

Ret/;-,' 



3 r m 

— I dcolm.[di(co - co) + 03(d) - co)]lms^ et (cd). 

K JO 



Differentiating the latter expression with respect to and tak- 
ing into account symmetries of integrands we obtain 



d 3 

-^-ImE = Im 

dA % 



day dd 3 
dco dco 



Im4 et (0) 



£0=0 



where ImE = Im-^-| ra= o. Note that the first derivative 

ImE = Im g° 1(3=0 vanishes. It means that (B-expansion 
of ImE Ret (o) « jlmT/'co 2 starts from the quadratic term, in 
compliance with the Fermi-liquid expression ( 1321 ). 

In order to find an fRG estimate for ImE", we need to es- 
tablish equations for a\ = lm^\m=o and a' 3 = Im^l^o. 
This is achieved by differentiating ( 11061 ) and ( 11071 ) 

da' U 2 

= -Ua[F 3 (0) + T lm/r 3 , (0) ) (152) 



U 2 . 



^ = Ua' 3 F 3 (0) + ^-MF:M)). 



(153) 



Using ( 11231 ) and Ims Ret (0) = 2/T?, we obtain the following 
set of equations 



dlmE" 
dA 
dd[ 
dA 
da' 3 
"dX 



= -(a[+a' 3 ) 



_ f 4U 4U 2 



%T 2 



4U 4U 2 

-a 3 



The solution of the last equations is given by 



4U -4£L 
1 e ^ 



7rr A 
4U 



^7 



(154) 
(155) 
(156) 

(157) 
(158) 
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and therefore 



known from Bethe-Ansatz and are 50 



1 — cosh- 



AU 



(159) 



Integrating the equation for ImE", we obtain at A = 



ImE" = - 

r 



, nT . , AU 

1 sinh — 

AU kT 



(160) 



At large U T an exponentially large scale emerges in a' 3 ss 

AU i AU 

— je*r and in ImE" « -w eW ■ We a l so note tnat fl 3 deter- 
mines the width B = —U /a 3 given in Eq. ( 11281 l. 

In a similar way we can establish the Fermi-liquid coeffi- 
cient corresponding to small-temperature expansion. Let us 
set co = in Eq. (I150l l and rescale the integration variable 
CO =xT. Then 



dA 



ImEHO) 
dco 



(161) 



n J sinh Co /T 
3T f dx 



Im [fl] (cd) + 03(d))] Ims^d)) 



Im [oi (xT) + a 3 {xT)} lms^(xT). 



n J sinhx 

Differentiating it twice with respect to T we obtain 



d d 2 Eg et (0) 
dA Im dT 2 



(162) 



T=0 



2xdx 
sinhx 



Im 



dai daj, 
dco dco 



Im4 et (0), 



co=0 



where f"^ 3^ = n 2 . Comparing this expression with ( 1154b 
we establish [cf. Eq. ( f32b l 

a 2 E* et 



Im- 



«9r 2 



= 7T 2 ImE". 



(163) 



m,T=Q 



Calculation of the Fermi-liquid coefficient corresponding to 
small-bias expansion is presented in the AppendixlDl 



Zs = exp(|^) 



U Jo ^"^"'^ 



cos(7Dc/2) 



(164a) 



Xc 



= «P("f?) 



r f , f kF 2\ COSh(7Dc/2) 

-y o Ckexp(-^ x )___. 



(164b) 



For U > TzT the reduced susceptibilities acquire their asymp- 
totic form where^S 



Xs i Xc — Xs 



«p(f [r-F])- ^* r - 

(165) 



This exponential behaviour governs the width of the Kondo 
resonance, which can be seen as follows. At eV s = (= jJ,), 
T =0, B = 0, V = the retarded Green function can be ap- 
proximated for co close to by 



Gf l (co) 



1 



(166) 



m*0) + ir/2' 

compare (01 and d32l . the effective mass m* being defined as 
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Hence 
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so that the full width at half maximum of the peak within this 
approximation is given by 
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From ( l32t and ( 11651 ) we deduce 
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VIII. RESULTS OF THE FREQUENCY DEPENDENT FRG and thus 



In this section we discuss numerical results obtained from 
the frequency dependent Keldysh fRG and compare them to 
those found by other methods. The underlying approximation 
scheme is the one described in section [VTl in particular, only 
the static part of the self-energy is fed back into the flow of E' 1 
and 7j , as explained in section lVI.DI In the figures, results of 
this method are labelled by "fRG". 

We start with the particle-hole symmetric model in ther- 
mal equilibrium at B = 0. The low energy expansion of the 
retarded self-energy is then given by Eq. (|32l . The reduced 
susceptibilities serving as coefficients of this expansion are 
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Figure [2] shows a comparison of the effective mass com- 
puted by the fRG to that obtained from second order pertur- 
bation theory and to the Bethe-Ansatz result. Since second 
order perturbation theory will serve us several times as a ref- 
erence solution we briefly sketch how it is determined. When 
evaluating the frequency-dependent second order diagram for 
the self-energy we use propagators dressed with the restricted 
self-consistent Hartree-Fock self-energy as internal propaga- 
tors. The self-consistent Hartree-Fock solution is obtained 
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FIG. 2: The effective mass m* = 1 — Hi 



as function of U for 



-a | £o=o ' 

V g = 0, T = 0, B = 0, V = 0. The exponential behaviour exhib- 
ited by the Bethe-Ansatz result obtained from Eqs. d!70ab and {T64J 
at large U /T is neither reproduced by the fRG nor by second or- 
der perturbation theory. While the effective mass computed from the 
frequency dependent fRG with static self-energy feedback (labelled 
"fRG" and used as default throughout this section) is comparable 
to the one found in second order perturbation theory, the values ob- 
tained from the fRG with full self-energy feedback are clearly worse 
(compare section PVI.Db . 



from the self-consistency equation 



-Ret 



U{n<f) = -—U jdaG^a)) 



= ^/da,/ eff (o)) (Q) _ % _^ et)2+r2/4 . 



(172a) 
(172b) 



This equation has a unique solution for small U, but there 
are three solutions at larger values of U (e.g. for U > nT/2 
for T = 0, B = 0). Two of them feature a local moment, 
E^ et 7^ T.f et , even for B — 0. 43 The notion "restricted Hartree- 
Fock" refers to the third solution which is energetically unfa- 
vorable as compared to the other two ones but is nevertheless 
preferable because of preserving spin symmetry, E^ et = S^ et . 
In the particle-hole symmetric case, for instance, the restricted 
Hartree-Fock self-energy E^ et = T,^ et = [7/2 aligns the single- 
particle levels with the chemical potential. We use this pertur- 
bation theory only for B = 0. For finite magnetic field it is not 
clear which solution of the self-consistency equation ought to 
be used to dress the single-particle propagators. 

In Fig. |2] we see that the effective mass computed from the 
fRG does not increase exponentially for large U, but it is com- 
parable to the one computed with second order perturbation 
theory. As a consequence the width of the Kondo resonance 
of the spectral function resulting from the fRG approximation 
does not shrink exponentially for increasing U. Figure [2] also 
shows the effective mass obtained from the fRG with full fre- 
quency dependent back-coupling of the self-energy into the 



FIG. 3: The effective mass as function of U for Vg = 0, T = 0, 
6 = 0, V = 0. The figure compares the Bethe-Ansatz and fRG 
data of Fig. [2] with results of the Matsubara fRG in "approximation 
1" of Ref. |57|, the Matsubara fRG with partial bosonization in the 
spin-rotational invariant version of Ref. 53 (data extracted from the 
preprint arXiv:081 1.2809 vl), and with the NRG data used for com- 
parison in Refs. 53.57. 



flow. As mentioned in section [VLD] this feedback leads to an 
artificial broadening of spectral features which becomes ap- 
parent here in the too low values for the effective mass. In 
Fig. [3] a comparison with the results of recent Matsubara fRG 
implementations is given. 

From the expansion ( l32l and from ( 1165b and (1170b it also 
follows that 



ir 
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U>7tT, (173) 
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increases exponentially with U in the Kondo regime. Figure|4] 
presents the respective data obtained by fRG and second or- 
der perturbation theory in comparison to the Bethe-Ansatz re- 
sult. In contrast to second order perturbation theory the fRG 
produces an exponential scale as has been proven already in 
Eq. ( 1160b . Using Eqs. ( 1169b and (1173b we can relate this ex- 
ponential behaviour to a Kondo scale 



\T 1 A UY ( 2 U 



(174) 



The identical exponential behaviour has been found to govern 
the pinning of the spectral weight to the chemical potential 
in the Matsubara fRG with frequency independent truncation 
schem e 54 ! 55 . With 2/n ~ 0.64 the prefactor is slightly smaller 
than the correct one in Eq. ( 11711 ). 7t/4 ~ 0.79. Figure [4] also 
shows that the fRG with feedback of the full frequency de- 
pendent self-energy into the flow does not produce an expo- 
nential scale in the second derivative of the self-energy. This 




FIG. 4: The second derivative of the self-energy as function of U 
for V g = 0, T = 0, B = 0, V = 0. The fRG produces an exponential 
increase with the exponent 41/ / nT, whereas the exact exponent is 
TtU /2T. Second order perturbation theory and the fRG with the full 
frequency dependent self-energy feedback (compare section IVI.Dt 
do not yield an exponential increase. 



illustrates again that the static feedback scheme which we use 
as standard is advantageous, compare section fVI.DI 

The fact that the fRG produces an exponential behaviour 
in ImE" but not in the effective mass constitutes an incon- 
sistency which manifests that the method captures Kondo 
physics only partly. The exact result satisfies 
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" Xc^ 



> 



(2s Xc) 



~2 
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(175) 
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where both sides asymptotically approach j£ s 2 /4. Figure|5]il- 
lustrates that the fRG result violates this inequality for U > 
4r. Therefore we restrict the following discussions to U < 4T. 

Figures [6HH] show spectral functions in the particle-hole 
symmetric case at B = 0, T = 0, in equilibrium. Results of the 
fRG and of second order perturbation theory are compared to 
essentially exact NRG data 77 . For U < 2Y all three methods 
yield nearly identical results. For U — 2F the shape of the 
peak differs already significantly from the Lorentzian form 
produced by static approximations such as restricted Hartree- 
Fock. At this interaction strength the difference occurs al- 
most exclusively due to the second order self-energy diagrams 
which are also captured exactly by the fRG. For U = 3r,4r 
the resonance peaks produced by fRG and second order per- 
turbation theory are too broad, see insets of Figs. [7] and [8] 
which is a consequence of the effective mass being too small 
as discussed above. The overall shape of the spectral func- 
tion is reproduced better by second order perturbation theory 
than by fRG The fRG has a tendency to shift spectral weight 
too far away from the central peak; in Fig. [8] the side peaks 
computed by fRG are situated at \co\ ~ 3r whereas they are 



FIG. 5: Comparison of the square of the effective mass and of the 
second derivative of the self-energy as function of U for Vg = 0, T = 
0, B = 0, V — 0. In the exact solution both have identical large (J 
asymptotics, m* 2 being always the larger one of both. In the fRG 
data the curves cross at U slightly above 4T. 
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FIG. 6: Spectral function for V g = 0, T = 0, B = 0, V = and U = 2T. 
The results of fRG and second order perturbation theory agree very 
well with NRG data. The peak shape differs already significantly 
from the Lorentzian form produced by static approximations such as 
restricted Hartree-Fock. 



expected to be at | £»| ~ U/2 = 2F. For the computation of 
equilibrium spectra it is hence preferable to resort to the fre- 
quency dependent Matsubara fRG which has been found to 
be clearly superior in this respect to second order perturbation 
theory at U = 2.5TM 

The linear conductance as function of gate voltage for dif- 
ferent magnetic fields at T = and U = 4T is shown in Fig. [9] 
The fRG data agree very well with NRG results; the frequency 
dependent Keldysh fRG obviously maintains the high qual- 
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FIG. 7: Spectral function for Vg = 0, T = 0, B = 0, V = and U = 3F 
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FIG. 9: Linear conductance as function of gate voltage for T = 0, 
U = 4r and different magnetic fields. 



0.5 




1.5 



eVJT 



FIG. 8: Spectral function for Vg = 0, T = 0, B = 0, V = and V = 
4r. The resonance at ft) = produced by fRG and by second order 
perturbation theory is too broad. The side peaks resulting from the 
fRG calculation are situated at too large |ft)|. 



ity in describing Gn n (V g ,B;T = 0) which has already been 
achieved by the fRG with frequency independent truncation 
schem o 54 ! 55 . 

Figures [10] and [TT] present the linear conductance as func- 
tion of gate voltage for different temperatures at B = and 
U = 2r,3r. In contrast to second order perturbation the- 
ory the fRG reproduces the NRG results for the width of the 
plateau at T = and the position of the maxima for T > 
quite accurately. The decrease of the conductance at V g = 
with increasing temperature is captured not completely cor- 
rect by the fRG but distinctively better than by second order 
perturbation theory. The fRG provides an altogether accept- 
able description of G\i n (V g ,T) for U < 3. This is an impor- 
tant improvement compared to previous fRG approaches. The 
static fRG 54 55 is in principle unable to reproduce the mini- 
mum of Gii n at Vg = 0. The reason is that the linear conduc- 



FIG. 10: Linear conductance as function of gate voltage for B - 
U = 2T and different temperatures. 



tance is given by 



Clin = e 



,r L r R 
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dm 

dco 
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(176) 



as follows from the current formula (l22l . In any static approx- 
imation at B = the spectral function p{co) is represented by 
a Lorentzian peak with fixed width and height, centered at the 
renormalized level position. Obviously Eq. (1176b can produce 
a single maximum of Gi; n (Vg) when the renormalized level is 
aligned with the chemical potential, but no local minimum. 
For the frequency dependent Matsubara fRG on the other 
hand, problems with the analytic continuation from imaginary 
to real frequencies obstructed the computation of the linear 
conductance. 57 Although it has been later found possible to 
circumvent this obstacle, the results for Gi; n (Vg) are less ac- 
curate than those obtained from the Keldysh fRG 78 A static 
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FIG. 1 1 : Linear conductance as function of gate voltage for B - 
U = 3r and different temperatures. 



fRG scheme based on a real frequency cut-off in Keldysh 
formalism has also been found to reproduce qualitatively the 
shape ofGiin(yg,r)^In view of the argument presented after 
Eq. ( 1176b this has to be the consequence of a renormalization 
of the level broadening. The latter however can be achieved 
only in a dynamic approximation as can be seen from the fluc- 
tuation dissipation theorem d94l i that manifestly requires a fre- 
quency dependent Keldysh component of the self-energy if 
E Ret has nonvanishing imaginary part. Therefore we suspect 
that the corresponding result of the static flow scheme used in 
Ref. [39] is an artifact connected to the violation of causality as 
consequence of the real frequency cut-off. 

The results of our Keldysh fRG for the current as function 
of bias voltage at V g = 0, B = 0, T = have been compared in 
Ref.[l3|to data obtained by a TD-DMRG treatment. Excellent 
agreement between both methods has been found for U < 4T. 
A more sensitive quantity for comparisons is however the dif- 
ferential conductance. Unfortunately, the TD-DMRG conduc- 
tance data existing by now are not enough for a meaningful 
comparison. In Fig. [12] we compare the differential conduc- 
tance as computed by fRG to the fourth order perturbation 
theory results of Ref. 59. The agreement for small and large 
voltages is rather good. A discrepancy is found at intermedi- 
ate V where the fRG does not reproduce the anomalous peak 
found by Ref.|59jfor sufficiently large U. Probably the anoma- 
lous peak is an artifact of fourth order perturbation theory; 
this conclusion is supported by recent nonequilibrium QMC 
results for the SIAM.— 

The data for Gdis{V) at different magnetic fields presented 
in Figs. |T3l and [T4l could in principle be compared to results 
of the recently introduced scattering states NRG 12 . The data 
of the latter approach is however is still too noisy to allow 
for definite conclusions. Finally, Figs. Q3] and [l6]present the 
differential conductance for different temperatures in compar- 
ison to second order perturbation theory. 

A detailed comparison of nonequilibrium results of the fRG 
with that of recent numerical methods is given in Ref. |6U 



FIG. 12: Differential conductance as function of bias voltage for 
V g = 0, B = and different values of U. The temperature is zero 
except for the two fRG curves with U = F,2T where T = 0.02r; 
the corresponding two data sets for the current at T = have not 
been sufficiently smooth to allow for numerical differentiation. The 
fourth order perturbation theory results have been extracted from the 
preprint arXiv:cond-mat/021 1616vl to Ref. |59|. 
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FIG. 13: FRG results for the differential conductance as function of 
bias voltage for V g = 0, T = 0, U = 3T and different values of the 
magnetic field. 



IX. CONCLUSION 

In this paper we studied a frequency dependent fRG ap- 
proximation for the Anderson impurity model. In equilibrium 
this model provides the possibility to compare our results to 
reliable data from Bethe-Ansatz or the NRG. Existing studies 
of the model within Matsubara fRG give indication on how 
to approach the problem within Keldysh formalism. The non- 
equilibrium properties of the model are subject of numerous 
investigations in recent time; however, a consistent picture did 
not yet emerge. 
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FIG. 14: FRG results for the differential conductance as function of 
bias voltage for V g = 0, T — 0, U = 4T and different values of the 
magnetic field. 



FIG. 16: FRG results for the differential conductance as function 
of bias voltage for V g = 0, B = 0, U = 4F and different values of 
temperature, in comparison with second order perturbation theory. 
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FIG. 15: FRG results for the differential conductance as function 
of bias voltage for Vg = 0, B = 0, U = 3T and different values of 
temperature, in comparison with second order perturbation theory. 



In order to preserve the Fermi liquid properties of the 
model we have chosen as flow parameter the level broaden- 
ing T/2 — A. In section[V]we studied the basic approximation 
scheme with static vertex renormalization and reproduced the 
results of the static Matsubara fRG. For the more advanced 
dynamic second order truncation scheme it was necessary to 
parameterize the frequency dependence of the vertex in order 
to obtain tractable equations. The self-energy which we feed 
back into the flow is static, and does not contain interaction 
induced quasiparticle decay rates. Nevertheless, it turned out 
that the results can be trusted up to values of the Coulomb in- 
teraction U < 3r = 6A, which are typical parameter regimes 
also used in diverse recent numerical studies of the model, see, 
e.g., Ref. [6l|. Furthermore, for technical reasons we restricted 
our study to 7L = Tr and Tl = Tr. 



At the particle-hole symmetric point the description of the 
spectral function obtained by our fRG is comparable to sec- 
ond order perturbation theory at weak interactions, while get- 
ting worse for increasing values of U. But in contrast to per- 
turbation theory the fRG generates exponential behaviour in 
certain quantities. The exponent of the Kondo temperature 
extracted from the second order derivative of the self-energy 
is identical to that found from the pinning mechanism in the 
static fRG, which is not far from the exact value. The fact that 
the effective mass derived from the fRG does not exhibit an 
exponential scale is related to the missing renormalization of 
the static part of the two-particle vertex in the present second- 
order truncation scheme. An extension which takes into ac- 
count contributions from the flow of the three-particle vertex 
could probably cure this deficiency. In that case we expect that 
also the feedback of the full frequency dependent self-energy 
into the flow will be favorable compared to the static feedback 
used here. 

A big success of the method is the very good description 
of the linear conductance as function of gate voltage, temper- 
ature and magnetic field. At zero temperature excellent con- 
ductance data are obtained even for U — 4T while for finite 
T acceptable quality is achieved for U <3T with much better 
results than in perturbation theory. In this respect the method 
is a utile complement to existing fRG approaches: the static 
fRG is applicable only at T — 0, and the frequency dependent 
Matsubara fRG does not produce conductance data of compa- 
rable quality. 

Concerning non-equilibrium properties the results of the 
fRG compare rather good to those of other methods for U < 
3r as far as comparisons are possible. In particular, we ob- 
tained very good agreement for the nonlinear current I(V) 
with TD-DMRG data at the most critical point V g = T = B = 
0. Together with the fact that our equilibrium results agree 
well with NRG at finite magnetic field and finite temperature, 
we expect our results for the nonlinear conductance G(V) to 



23 



be reliable at finite T and B, although no firm benchmark has 
yet been obtained by other techniques. 

An issue to be addressed for future research is the prob- 
lem of current conservation by the fRG. For the finite bias 
data which we presented in section IVHll current conservation 
is given as a consequence of the complete symmetry estab- 
lished by eV g = (jUl + /Jr)/2, T l = T R , T L = T R . In more 
general cases the fRG is not expected to be a current con- 
serving method. Similar problems are known for other meth- 
ods, for example perturbation theory 58 . Only when reliable 
benchmarks become available, it can be concluded, whether 
the violation of current conservation generates serious errors 
in the nonlinear conductance G(V), and whether more elab- 
orate parametrizations within nonequilibrium fRG are neces- 
sary. 

The comparison of the frequency dependent Keldysh and 
Matsubara fRG reveals that in spite of being closely related, 
the two methods are not equivalent in equilibrium. The choice 
of an imaginary frequency cut-off as used in Matsubara fRG is 
not possible in the Keldysh fRG without destroying the Fermi 
liquid property at the particle-hole symmetric point. For the 
Keldysh fRG the static self-energy feedback into the flow is 
preferable to the feedback of the full one, as opposed to the 
Matsubara fRG And while the Matsubara fRG yields better 
results for the form of the spectral function at the particle- 
hole symmetric point, the Keldysh fRG produces better data 
for the linear conductance as function of gate voltage, mak- 
ing accessible finite temperatures where the Matsubara fRG is 
handicapped by the problem of analytic continuation. 

In total our investigation gave insight on how a frequency 
dependent approximation to the Keldysh fRG can be con- 
structed, and produced a flexible method to compute proper- 
ties of the Anderson impurity model for moderate interactions. 
A fundamental problem of applying the fRG in the considered 
truncation scheme to this model is that the justification for the 
performed approximations can only be given by perturbative 
arguments. This difficulty is rooted in the regular perturba- 
tive expansion of the model. This contrasts with the situation 
in models with low energy divergencies which typically allow 
to distinguish relevant and irrelevant contributions to the RG 
flow by power counting. It would be interesting to figure out 
whether perturbative RG techniques are still capable to cap- 
ture the exponentially small scale of the Kondo temperature in 
the effective mass for large Coulomb interactions. Within our 
proposed approximation scheme, this scale already appeared 
in the second order derivative of the self-energy, showing that 
our perturbative RG scheme has the potential to extract such 
scales. Therefore, it might appear possible that certain con- 
tributions from the 3 -particle vertex can help to reveal the 
same scale for the effective mass as well. We expect that the 
approximation can be systematically improved by including 
the influence of higher order vertices on the flow, if care is 
taken to formulate the approximation scheme in accordance 
with Ward identities. It might be also interesting to access the 
Kondo limit of the Anderson model using an approach based 
on Hubbard-Stratonovich fields decomposition as proposed in 
Ref. El in combination with nonequilibrium fRG 

In summary, the method proposed in the present paper 



seems to be sufficient to describe Coulomb interactions up to 
the order of the band width D of the local system. For a single- 
level quantum dot, the band width is given by D ~ T, whereas 
in quantum wires it is given by D ~ t, where t is the hop- 
ping matrix element. In both cases, the nonequilibrium fRG 
schemes of this paper and of Ref. [25] seem to be reliable for 
U <D, providing the hope for a unified approach describing 
the whole crossover from local to extended quantum systems. 
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Appendix A: Identification of independent components of <p^' x '' 



The numerous Keldysh and spin components of the func- 
tions <p?' x,d are not independent of each other. In this sec- 
tion we identify a set of independent components from which 
all the others can be determined. For that purpose we make 
extensive use of the relations for the vertex functions which 
originate from permutation of particles, complex conjuga- 
tion, causality and the KMS conditions and are described in 
Ref. |42[ As explained in that reference, the hybridization 
flow parameter preserves these relations during the truncated 
RG flow. This is as well the case under the additional ap- 
proximation induced by the replacement ( TTTT i since the static 
renormalized vertices appearing in the flow equation ( T8Tb have 
the same spin and Keldysh structure as the original bare ver- 
tex. The relations can be even used for the three individual 
functions <p?' x,d since the diagrammatic structure of the three 
corresponding flow equations is invariant under the manipu- 
lations done in the proof of those relations, cf. Ref. The 
only exception is that permutation of either the incoming or 
the outgoing indices maps <p x onto <p d and vice versa. 
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1. Analysis of the spin structure 

Due to the spin structure of v, (4>^' x ' d ), and (/P p ' xph ) certain 
spin components of <pf ,x ' vanish: Knowing that 

v o[o' 1 \a\<j2 = 
(^l)a[a^a 2 = }if <Ti ^ W 2 or c[ ^ g' 2i 



°2 I °2 



= 



(^ d ) CT ;^|<r lCT2 - if of + a' 2 =£ a, + ct 2 , 



(Ala) 
(Alb) 



we conclude from the flow equation d8TT i that 

(<Pl)a;(T|| CTl a 2 = 0, if CTi 7^ ^2 or a; ^ a' 2 , 
(<Px)a 1 '(T|| CTl a 2 =0. if(Ti + C72^ffi , + 0-2, 
(<Px)(T 1 '^| CTl (T 2 =0. if °1 + °2 £ Oi+oJ. 



(A2a) 
(A2b) 
(A2c) 



Concerning the remaining spin components, we can restrict 
our study to 
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because the other ones can be derived from these by permuta- 
tions of particle indices, 
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2. Analysis of the Keldysh structure 

In the following it is convenient to describe the Keldysh 
components of two-particle functions in terms of block matri- 
ces, in which they are arranged according to the table 



/(11|11) (11|21) 
(21|H) (21|21) 



(12|11) (12|21) 
\ (22|11) (22|21) 



(11|12) (11|22)\ 
(21 1 12) (21 122) 1 



(12| 12) (12|22) 
(22j 12) (22|22)/ 



(A5) 



The indices tt' 2l U2 of an index tuple (a[(X^\(XilXi) indicate 
which of the blocks is to be chosen, while a[ , (X\ defines the 
position inside a block. The bare interaction vertex (l40b for 
example is described by the block matrix 



(A6) 





(0 1 


1 0\ 


_a[a! 2 \a\Ui 1_ 


1 


1 


V ^tT||<Tl<T 2 ~~ 2 V<T l' <T 2l CTl °' 2 


1 


1 




Ko i 


1 0/ 



From the flow equation ( 18 lab it follows that <pf (n) has the 
Keldysh structure 



<Pa 



(c 


d 


d 


C \ 


a 


b 


b 


a 


a 


b 


b 


a 


\c 


d 


d 


c) 



(A7) 



which can be seen as follows: the initial value <p? =00 = is 
consistent with the structure (IA71 >. Since q>? (TV) has this struc- 
ture, so does [v+ <P?(n) +'&\ +*a1- P rom Eq. d81a| ) it can 
then be derived that d<pj /dA has the structure jA7| > as well. 
Therefore this structure is conserved during the flow. 

Exploiting the causality relation (<p£ ) 22 l 22 = and the trans- 
formation properties with respect to complex conjugation we 
find 



(<plW<r(n) 






a x 


p * 

a\ 





a\ 


b\ 


b\ 


a x 


a X 




b\ 


a x 





P * 

a x 


P * 

"A 






(n), (A8) 



where (a P jJ a a(n) is a complex valued function while 
{b\)aa(R) is purely imaginary. As a consequence of causal- 
ity, (aj^cjo^n) is analytic in the upper half plane of IX The 
relations for exchange of particle indices yield 

KW( n ) = ( a l 
K W(n) = (blha 



4W(n) 



(A9a) 
(A9b) 



In case of thermal equilibrium at temperature T and chemical 
potential = the only nontrivial component of the general- 
ized fluctuation dissipation theorem of Ref.l42l reads 



(tfW(n) 



i2coth( — ) Im(4)^(n). 



(A10) 



The same reasoning as used for {<p^) a o\oo(Vl) shows that 
the Keldysh structure of (9? )aa\aaQ0 i s 



(<Px)aa\oo(X) 








fl A 




a x 


b \ 





«r 


«r 





K 


"A 




"A 








(X), (All) 



with a complex valued function (af ) CT ct(^) and a purely imag- 
inary function (M)ctct(.X'). ( a ^)aa(^) is analytic in the lower 
half plane of X. The behaviour of a\ and Z?^ under exchange 
of particle indices is 

(4)^W = (4WK)*. ( A12a ) 

(^)^(X) = (^W(-X). (A12b) 
In case of thermal equilibrium the KMS conditions are 

(b\) aW (X) = -i2coth (£\ lm(a\)^(X). (A13) 
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The Keldysh structure of {(pf)a'S\&s(^) nas the form 



Oi)<t(t|<tct(A) 








(4%a 


(b{)oo 







(b{)oo 











(4W 












(A), 
(A14) 



where (a^) a -^(A) and (^^-^(A) are complex valued func- 
tions which satisfy 

(4W(A) = (aJW-A)*, (A15a) 
(ijt)«rtf(A) = -(fti)« rt f(-A)* J (A15b) 
(*l)aa(A) = (*lW(-A). (A15c) 

The function (a d ^-^(A) is analytic in the upper half plane of 
A. In case of thermal equilibrium the KMS conditions demand 



where (4 ) <t<t (A) is a complex valued and (b^ ) aa (A) a purely 
imaginary function; they satisfy 



(fljt)aa(A) = (fli)aa(-A)* 
(&£)aa(A) = (fti)aa(-A). 



(A18a) 
(A18b) 



(4)cc(A) is analytic in the upper half plane of A. In case of 
thermal equilibrium the KMS conditions are 



(A19) 



(b%)aa{A) = i2coth ( — ) Im(4) ff(T (A). 



Re(^)^(A)=tanh( — ] Re 



(4W(A)-(4)„ a (A) 



(A16a) 



Im(^) a -(A)=coth( — ] Im 



A 
2T 



(4)^(A) + (4)- ff ( A ) 



(A16b) 



Finally the Keldysh structure of ((p^) a a\aa(^) i s 



(<Px)oo\oo(£) = 


^ 4 


4* 4 ^ 
4 4* 


(A), (A17) 

I 


4* 4 
, wl „d* 


^ 

„d n 


V A A 





In summary we conclude that all spin and Keldysh com- 
ponents of the functions (pP x d can be determined from the 
selection 



(afV (^' x ' d ) n , (4) ;t , (alU, (4) 



(A20) 



The precise form of the flow equations shows, that (4)lt can 
be expressed through (4)t4-> see ^1- 6^ below. 



Appendix B: Flow equation for the self-energy 



In our approximation the two-particle vertex acquires a special structure described in section [VI. Bl and Appendix [A] When 
we insert this structure into Eq. (f2Tb we obtain as flow equation for the retarded component of the self-energy 



dl^ a (co) _ i 
dA ~ ~2n 



U 



dco' 



(ft* Wfl>' - 0)) + (4)^(0)] S§ et V) - [(4)^(0) - co') (4) CT<T (0)] S* et V) 
(4)^(0)' + a)) + (4)^(0)1 (a)') + (4) CTCT (0)Sf x {co') 



-Z + (fl^)aa(w' + (0) + (4) CT <f(<»' - O) + (4 



4 A (a)')- (4)<x<x(a>-a>')-(4)aa(0) 



(Bl) 



As a consequence of Eq. d36l > S Ret ((o') [S Av ((o')] is analytic in the upper [lower] half plane of ft)' and vanishes as ft)' 2 for 
I ft)' I — > °°. Hence, 



Jdco'S ReuAv (co')=0, 



(B2) 
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and Eq. ( IBU is reduced to 



dE ^ (6)) = jda>'\ {b\U{co> co)S™ x (co') (b{) aa {co ft/)S* et V) + {b\)^{co' + a>)S$* « 



2 + {"DM®' + ( °) + (4W(V - £») + (flj W(0) 



(B3) 



Since the relation E^ v (ft)) = E R (ft))* is maintained during the flow, it is not necessary to compute the flow of E^ v (ft)) 
separately. Following similar steps as for E Ret we acquire the flow equation for E , 



dE£ A (ft)) i 



dA 



2% 



dco' I 2 Re 



| + (aj[ Wa>' + ft))* + (4 W(»' - «) ) 4 et " 



2Re 



(4) CTCT (ft)-o)')S* et V) 

+ [{b\) a a{C0' + ft)) + (^)^(ft)' - ft))] SfV) - (ii) ffff (fl) - (0')Sl X {<0')\. (B4) 



Appendix C: Flow equation for the two-particle vertex 



Due to the structure of the two-particle vertex function described in section IVI.BI and Appendix [A] its flow is determined 
completely by the flow of the components ( IA20b . The flow equations of these components can be formulated as a closed set, if 
the other components are eliminated by use of the relations found in Appendix[A] Setting up these flow equations is cumbersome 
but straightforward. This section is devoted to a brief sketch of some simplifications which can be made during the calculation 
and indicating the result. 

The flow equations for {d^aa = (<Pf ancl (b\)cro = f°ll° w fromEq. (18 lab . Due toEqs. {79]), d82|, and ( lA2ab 

the implicit summation over spin indices in Eq. (18 lab is reduced to the two contributions ( <T3 CJ4 1 CT4) = (cra\<ja), ((7(7 1 (7(7). 
The summation over Keldysh indices is also largely reduced: first, because of the vanishing components in Eq. dA8l ), and, 
second, because of 



(/PP)" 1 " 2 ^ =0, if «i = «i = 1 or a 2 = c£ = 1, 



(/ P P) 12| 21 = (/f) 



p P) 21|12 = Q _ 



(Cla) 
(Clb) 



Equation dCTab is a consequence of G 1 ! 1 = 0, S 1 ! 1 = 0, while Eq. ( fClbt follows from G Ret (ft)), S Ret (ft)) [G Av (ft)), S Ay (ft))] being 
analytic in the upper [lower] half plane of ft). Making further use of 



2(g)) = _(_l)«l+«2+»l + a 2(/PP-)«I«2|ai<^( (0 -)* 



(C2) 



for the remaining components, which follows form G Ret (ft))* = G^ v (ft)), S Ret (fi))* = S Av (©) and G^(ft)) 
-Sf(fi>), we find 



-G K (ft)),S K (ft)) 



d(g|W(n) 

dA 



u+m+ 



• + KW(n) 



(/D 2 i 12 (n) + (/f) 2 i 21 (n) 



(C3a) 



W(n) 



dA 



= 2ilm 



KW(n) 



(/r)S n (n) + i(/f) 2 CT f 2 (n) 



+ 



+ (4 WW 



W(n) (TC 12 (n)+(/f)X(n 



,22121 , 



(C3b) 



The flow equations for {a\ ) a o{X) and (bV)aa(X) can be derived in an analogous way. Instead of Eqs. ( ICU and (IC2b we use 



and 



(/P h )«i«2l«!«2 = 0, if «i = a( = 1 or 02 = <4 = 1, 
^11,22 = (/ ph )22 |n = 05 



(C4a) 
(C4b) 

(C5) 
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and find 



dA 
dA 



u + u\ + u A 



+(4W( X ) 



« 22 (*)+(OS 12 (*) 



(C6a) 



= 2ilm 



£/ + £/? + C/ d 2 

2 — 



+ 



u + ul + u* 



(C6b) 



The flow equations for (a^o-o^A) and (b^o^A) are 



d(4)^(A) _ 


"C/H 




+ ui 


dA 




2 




dA 




2 


+ ui 








+ u* 






2 






c/- 




+ ui 






2 














2 



w d 



X a 



w, d 

(4W(A) — 



+ (4W( A ) 

+(4W( A ) 
+(4W( A ) 

+ (4W(A) 
+ (^)aa(A) S 



E 



(a*)„(A)--i* 



(/f) 2 . 2 l 21 (A) + (/f)i 2 ' 22 (A)' 



(C7a) 



(4)^(A)- 



/f)S 21 (A) + (/rff(A) + (/f)S 22 (A) 



'(/f (A) + (Jf )aci 12 (A) + (/f ) 22 f (A) 



(*1W(A) [(/f )S 21 (A) + (/, ph )# 2 (A) 
(fti)aa(A) [(/f C(A) + (jf ) 2 f 2 (A) 



(fciWA) f(/f) 2 4 22 (A) + (/f) 2 4 12 (A) 



(^)^(A) (/f)S 21 (A) + (/f%r(A) 



phs!2|22. 



Since ( 

^3 )(jff and )o 7 (j have the same (zero) initial value at X — °° we conclude from Ec[. dC7al ) that 

(4)aa(A) = (4W(A). 
Finally, the flow equations for (4)<t<t(A) and (#[) a(T (A) turn out to be 



(C7b) 



(C8) 



d(4)qg(A) 

dA 



— \ — ^+(4wa) 



(^ ph ) 2 4 2i (A)+(/r)§ 22 (A) 



(fll)aa(A)- 



A cr 



d(4MA) 
dA 



2ilm 



u + ul + ui 



(7f)^ 1 (A) + (7f),r(A) 

2 



(C9a) 



(A) 



(^ ph )S 2i (A)+^(/rff(A) 



(4w(A)-^a 



2 r 



(7f)^ 21 (A) + i(/f) 22 l 22 (A) 



+ (4W(A)* 



(*SW(A) (/f)S 22 (A) + (/f)^(A) 



phs22|12. 



w 



X a 



(bi) aa (A) [(/f ) 21 f (A) + (/f) 22 f(A) 



(C9b) 



Appendix D: Nonequilibrium Fermi-liquid coefficient We remind that these conditions are important to ensure the 

In this evaluation we make the same assumptions as in the 
section |Vni i.e. T L = F R and ji L = -jx R = sV,B = eV g = 0. 
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I / eV\ 1 / eV\ 
—hp ( © — — ] + —hp ( © + — ) , (Dl) 



real value of the static component of the vertex as well as the 
current conservation. 

Let us define a nonequilibrium distribution function 

h n F eq {(o) = a L h F)L ((o) + a R h F . R ((o) 
eV\ 1 , 

Y) + 2 h \- 2 J 

where ai + a R = 1 . 

The nonequilibrium equation for the self-energy (IBll l in the 
particle-hole symmetric case reads 

^E^(ft)) = ~JdG) [- 7o /4 e >)Im4 et (ft)) (D2) 

+h F e9 (&)a 1 (& + ft)) Im4 et (ftj) +bi(8) + ©)4 V (©) 
+h F eq (a>)a 3 (G3- ft))Im4 et (d)) +fe 3 (ft) - ft>)4 et (<») 

where the vertex functions a and b components are no longer 
related to each other by the KMS equalities, and where bia is 
defined on the analogy of Eq. ( 1146b . 

In the following we will need an expression for a product 
of the two distribution functions ( 1DU evaluated at different 
frequency arguments ft>i and ©2. At small eV such a product 
can be represented as 

[a L h FL (coi) + a R h Ffi (coi)} [a L h FiL {ayi) + a R h FR (a>2)} 
= a L h F , L (eoi)h FL {(»2) + a R h F . R ((Oi)h Ffi {o)2) 
-a L a R [h FtL (coi) ~ h FR {coi)} [hp >L {oh) -h F:R {a>i)) 
ss a L h FtL (a>i)h FL (o}2) + a^p.nico^hp^coj) 
-4a L a R (ey) 2 5(ft)i)5(ft>2). (D3) 

Let us consequently analyze voltage-induced corrections to 
the KMS relations occurring in each channel. 



1. Particle-particle channel 

In order to evaluate (11151 l we set a>\ = y + ft)', ©2 = § - ©' 
and use for the PP channel the identity 

h F , n (©1 )h F , n (©2) = -1 +h B , r] (©1 + ©2) [h F , n (©1 ) + h FjV (©2)] • 

(D4) 

where /! Bj) (n) = h B (H — rjeV), rj = ± for tj = L/R, and 
h B {co) = coth = sign(©). It leads to 

E <hh F>n (a)i)h Ftn (a>2,) = 

T)=± 

= -i+Jl^(n) [/zr^O+C'M 

+x X! 7 ? /! B,i?( n ){° : i[ /j f,i( G) i) + /j / ; 'x(aJ2)] 

^ 7]=± 

- ccr [h F>R ( ©1 ) + h FlR (©2 )] } , 



(D5) 



where 



£!1 ^(n)_ 2ey5(n)j 

77 = ± z 



(D6) 



as well as 

cc L [h F> L(oh) +h FiL (G>2)] - a R [h FR (a>i)+h Ffi (a)2)} 
«-eV[S(©i) + S(©2)] (D7) 

for (Xl = CC R = \. Therefore, the second sum in ( lD5t approxi- 
mately equals k, 4(eV) 2 <5(©')<5(n). 
We can represent 



// 1 (n) = //f(n)+//f / (n), 



(D8) 



where 



fff(n) = 
Hf(n) « »K+*?)*(n), 



E W n ) 



[F!(n)-F*(n)], (D9) 

(D10) 



and 



aLa«K) 2 ^-[lm^ et (0)] 2 , (Dll) 



-44. 



Therefore the solution of Eq. (Ill lb has the form 

b 1 (n) = b[(u)+i(b[ Ia +b[ Ib )8(u), 

where 

1 



(D12) 



(D13) 



E W n ) 

TJ=± 



[fli(n)-aKn)], (Di4) 



and fe| /cl and correspond to the inhomogeneity terms 
{eV) 2 in Eqs. dP3t and dD5l ), respectively. They can be found 
from the equation 



db 



II,a/b 



dA 



= \a l (0)\ 2 s a 1 /b +2u 1 (0)F l (0)b"' a/b . (D15) 



2. Particle-hole channel 



In order to evaluate (11161 > we set ©1 = ©' — §, ©2 = ©' + § 
and use for the PH channel the identity 

/!f, 7 ,(©l)/!F^(© 2 ) = I - h B ((Oi - 0)2)[h Fll (C0i) - h Fri ((02)} . 

(D16) 

Then 



H 3 (X)=Hl(X)+Hl I (X), 



(D17) 



where 



and 



= -h B (X){F 3 (X)-Ff(X)}, (D18) 
# 3 7/ (X) = is 3 S(X), (D19) 



. 3 = -a L a fl ( e y) 2 A [lm^ et (0)] 2 • (D20) 

7T Q.A 
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The solution of Eq. ( II 12b has the form 

b 3 (X)=b' 3 (X) + ib I 3 I 8(X), 

where 

b I 3 (X) = -h B (X){a 3 (X)-a* 3 (X)}, 
and b^(X) obeys the equation 
db 11 

= |fl 3 (0) | 2 ^ 3 + 2a 3 (0)/^ (0)^ 7 . 

Analogously, the function H 2a (A) occurring 
Eqs. (1137b and ( 1138b can be represented as 



where 



H 2a {A)=H I 2a (A)+H I 2 ' a {A), 



H' 2a (A) = h B (A){F 2a (A)-F2 a (A)}, 
H» a (A) = is 2a 8(A), 



and 



sia = -a L a R {eV) 2 — [lmg£' 
% dk 



Ret/ 



(D21) 
(D22) 

(D23) 
in the 

(D24) 



(D25) 
(D26) 



(D27) 



Note that from ( 11241 ) it follows 

si = S3 = S2c = -(eV) 2 Im^(0). (D28) 
The solution of Eqs. ( 1137b and ( 1138b then takes the form 
^20/2a(A) = b' 20/2a (A) + ib' 2 ' 0/2a 8(A), (D29) 

where 

b' 20/2a (A) = h B (A){F 20/2a (A) - F* 0/2a (A)}, (D30) 

and 



db 11 

= -2 fl20 (0)F 2(7 (0)^ / (J 



(D31) 



d/j 



// 

2a 



dA 



= -|fl2o(0)| 2 S2tT-2a 2 o(0) J F 2ff (0)fe 2 / . (D32) 



Comparing the latter relations with Eqs. dD 151 ) and dD23l ) and 
using (ID28b we establish that 



b" 



2a - 



■ 1 , 3 , b&= 3 - 1 . (D33) 



3. Self-energy 



An equation for the imaginary part of the self-energy in nonequilibrium reads 



— / d(0 
% 



h n ; q {6i))lma l (a))lms R ?\a))--b[(a))lm4^6i)+^ 



+ _L [(*,«« + fcf )Ims Av (0) + (b'j - Olm, Ret (0)] , 



which simplifies due to the particle-hole symmetry to the form 



(D34) 



dA ImE ^ 



— / da>< 

% , 



■rt=± 



Imai(a))Ims§- et (<» 



[h n p q {&i) - h B (6))] Im[a 3 (ffl) + a 2a {ob)} Ims^Xeb) 



1 

1% 



-2b\ la ' + ^(b\ la +b 3 ') 



Ims Ket (0) 



(D35) 
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Expanding it in eV, we obtain 



d ImE Ret (0) = " 



dA 



eV\ 



2% \ 2 ) 
1 



dco8\co)lm\d l {co)+a 3 {co)]\ms^Hco)--{eVY / rfffl5'(ffl)Imai(ffl)Ini5 Ket (ffl) 

7T ' 



-ib'^-bf + ^b'^+b'!) 



Ims Ket (0) 



da\ daj, 
<?<» dco 



Ims Ket (0) + -(eVrim 



fi)=0 



<9gm 



<9(B 



Ims Ket (0) 



£0=0 



1 
1% 



-2b 



Ila ,IIb 



Ims Ret (0) 



r 



(D36) 



Comparing Eqs. (ID15b and ( ID23t with d 152b and (1153b we After comparison with ( 1 1541 ) it becomes clear that 
establish that 





-(eV) 2 Im 


da\ 
dco 


03=0 


(D37) 


fcf = 


-4b[ Ia , 






(D38) 


ft? = 


-(eV) 2 Im 


<3«3 
d CO 


0)=0 


(D39) 



These relations allow us to express ( lD36b as 
d 



ImE Ret (0) 



dA 



3 3(eV) 2 
27T 4 



1m 



(9© <?<» 



Im- 



«9 2 E 



<?V 2 



-ImE" 



<B,V=0 



Ims Ret (0) (D40) which agrees with (|32j and Ref. [ 



<a=0 



(D41) 
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